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Description 

[0001] A portion of the disclosure of this patent document contains material which is subject to copyright protection. 
The copyright owner has no objection to the facsimile reproduction by anyone of the patent document or the the patent 
s disclosure, as it appears in the Patent and Trademark Office patent file or records, but otherwise reserves all copyright 
rights whatsoever. 

[0002] This invention relates to color matching and more particularly to a process and related apparatus for matching 
the color displayed on a plurality of color display devices. 

10 2. Description of Related Art. 

[0003] The primary use of this invention is in the field of computer assisted color publishing systems particularly in 
the area of color matching. In such systems typically a color image is scanned using a scanning device which measures 
light intensity reflected or transmitted from the image to generate electronic signals representing three achromatic 

is images of the original, each of the three achromatic images resulting from detecting the light intensity from the original 
color image through a colored filter. The filters used are almost always Red, Green and Blue. The electronic signal 
from the scanner is converted to a digital signal in which light intensity levels are represented as numbers. Information 
identifying each set of numbers representing the image information obtained through each filter is also preserved. 
Thus, through these steps, the original colored image is converted to a plurality of image values, and for each picture 

20 element in the picture, there are three such values, a Red, a Green and a Blue. 

[0004] The image represented by the image values may be displayed in a CRT type monitor, or may be printed using 
a printing device able to accept electronic input. Often the image will be displayed at different locations and times using 
more than one display device, such as a plurality of CRT displays for observation and study by more than one operator. 
Similarly, hard copies may be desired in two or more different locations using two or more different printers. Even 

25 though the input values to the multiple CRT displays are the same, the displayed colors are vastly different, as anyone 
is well aware of who has ever observed the multiple television displays in a typical television sales store front. However, 
when one needs to make decisions about acceptability of color for a display, one needs to know with a great degree 
of confidence that the color, that one is observing and discussing with an operator who is observing the same image 
on a different CRT, is the same as the color observed by the operator. The same is true, if the displayed image is one 

30 created on a printer and compared with the same image created on a different printer. 

[0005] Colorimetry, which is the study of color based on both spectral distribution of the energy reflected ortransmitted 
from a test sample and the response of the human eye, as well as the spectral distribution of the illuminating source, 
provides a method to describe and measure color and enables one to determine when colors match. Through the use 
of CIE defined Tristimulus Values (which are the amounts of three primary lights which when added produce a visual, 

35 or colorimetric match with an original color), one may determine with reasonable certainty that if two colors have the 
same three CIE Tristimulus values that is, if the Red Tristimulus value of one color is the same as the Red Tristimulus 
value of the other, and so on for the Green and Blue Tristimulus values, then the appearance of the two colors will be 
the same to the average observer. CIE stands for the International Commission on Illumination. 
[0006] It appears, therefore, that the problem of matching color outputs of different displays is readily resolved by 

40 providing displays having the same CIE tristimulus values for the same image value inputs. 

[0007] While the solution in principle is simple, creating a conversion device which will make two or more different 
displays produce the same tristimulus output values for the same input image values is very difficult. Each of the display 
devices operates in its own, device dependent, color space where image values at its input are transformed into display 
image values. The difficulty lies primarily in matching the two transformations occurring within the two display devices 

45 for the image values appearing at the input of each, so that both display the same displayed colorimetric values for 
the same input image values. 

[0008] The prior art solutions to color matching fall into two fundamentally distinct approaches. The first is based on 
the decomposition of a color vector to a set of primaries, and known as the primary decomposition technique. Typical 
of this approach is the technique described in U.S. Patent application 07/494,463 filed March 16, 1990, assigned to 

so E.I. du Pont de Nemours and Company and corresponding to EP-A-0 446 914. According to this technique, the input 
image is decomposed into the linear combination of eight primaries (yellow, magenta, cyan, black, red, green, blue, 
and three-color black). Coefficients for a 4x8 matrix (the values of CYMK for each of the eight primaries) are adjusted 
such that the color of the eight primaries in the input from both systems to be matched are matched, and the new input 
values are found by matrix multiplication of an 8X1 matrix (coefficient of decomposition) with the 4X8 matrix. This 

55 technique suffers because of the non linearity in the additive properties of color dyes and because the "primaries" used 
are not mathematically independent. 

[0009] The more common solution, adopted by the printing industry is the grid sampling technique. This involves 
using a color transformation formula such that the error between the targeted color and the processed color is mini- 
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mized. The differences between the various methods in existence are found in the specific transformation formulas. 
U.S. Patent No. US-A-4,500,919 issued to Schreiber is a good example of the use of a transformation formula. The 
difficulty with this approach again lies in the non-linearity of the color addition process and the complexity of the color 
surfaces in a set of equations. 

5 [0010] A method corresponding to the precharacterizing part of claim 1 is disclosed in EP-A-0 475 554. This reference 
describes a method for matching the color displays of a proofer and a printer. The method uses a printer conversion 
table which converts between a first input digital representation CMY and a measurable color coordinate system XYZ. 
A second conversion table is provided for the proofer. The contents of the conversion tables are delivered through a 
printerrto-proofer transformer to a transformation table for each entry in the printer conversion table, the operation 

10 involves searching the proofer conversion table for entries whose XYZ value is near to the printer XYZ. The results 
are stored in a "table of close values". The transformation table may be made more accurate by producing a new 
version of the proofer conversion table by measuring colorimetric values of prints that were printed by using the known 
CMY data. 

[0011] It is an object of the present invention to provide a method and an apparatus for matching the color display 
is of at least a first display device and at least a second display device, providing a high degree of accuracy of color 
matching. 

[0012] This method of the present invention is defined by claim 1 and the apparatus is defined by claim 9. 

[0013] The apparatus may further comprise means for storing the generated LUT. A gamut mapping means may be 

included between the output of the second model and the means to compare for mapping output values from the 

20 second model into a color gamut commensurate to that of the first model. 

[0014] The apparatus may comprise hardware, or may be a computer programmed through software to perform all 
of the 'above operations, or may be a combination of dedicated hardware and computer implemented software. 
[0015] The invention can be more fully understood from the following description thereof in connection with the ac- 
companying Figures described as follows. 

25 [0016] Figure 1 is a schematic representation of an arrangement in which two matched images are displayed on two 
CRT displays. 

[0017] Figure 2 is a schematic representation showing the arrangement used in obtaining data for generating the 
LUTs used in modeling the two displays of Figure 1 . 

[0018] Figure 3 is a schematic representation of the process for obtaining the transform LUT used in the adaptor in 
30 the arrangement shown in Figure 1 . 

[0019] Figure 4 is a generic representation of a display model. 

[0020] Figure 5 is a schematic representation of the process used in deriving the correction factors shown in Figure 3. 
[0021] Figure 6 is a schematic representation of an arrangement in accordance with this invention in which two 
matched images are printed using two printing devices. 
35 [0022] Figure 7 is a schematic representation showing the arrangement used in obtaining data for generating the 
LUTs used in modeling the two printing devices of Figure 6. 

[0023] Figure 8 is a schematic representation of the process for obtaining the transform LUT used in the adaptor in 
the arrangement shown in Figure 6. 

[0024] Figure 9 is a schematic representation of an arrangement in accordance with this invention in which two 
40 matched images are printed using a CRT display and a printer having including an RGB to CMYK converter. 

[0025] Figure 10 is a schematic representation of an arrangement in accordance with this invention in which two 
matched images are printed using a CRT display and a color proofing device using a YMCK to RGB to drive the CRT 
display. 

[0026] Figure 11 is a generic representation of a display model for use in developing the transfer LUTs required for 

45 the embodiments shown in Figures 9 and 10. 

[0027] The invention will next be described in detail with reference to the drawings in which similar characters indicate 
similar elements in all figures of the drawings. During the following discussion, color values are quantized in 256 steps 
corresponding to an 8 bit system. Other quantizations are possible, and not intended to be excluded because of the 
use of the 256 steps in the following examples. 

so [0028] The invention comprises a method and apparatus for providing color visual matching of two representations 
of the same image when the image is displayed in two different displays, so that an observer will on visual examination 
of the displayed images observe two images having substantially the same color appearance. Using the method or 
apparatus of this invention causes each of two observers, whose visual response substantially conforms to the visual 
response of an average observer as defined by the CIE institute, upon observing any one of the color displays under 

55 the same or similar surrounding illumination and background, to receive substantially the same visual impression. 
[0029] Figure 1 shows in schematic representation a situation where a color image is displayed on two CRT type 
displays which receive RGB type inputs to display an image. The image may be stored and manipulated in digital 
format in a work station 1 0 which itself may comprise a scanner, an image processor, a display device, and have image 
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storage capabilities, or any portion of the above, or more. 

[0030] Because the required input to displays 12 and 14 in this example is a digital RGB type signal, the workstation 
outputs over line 16 a digital RGB signal representing a color image. The digital RGB signal comprises a set of color 
values, one set of three values (R, G, B) for each picture element of the image to be displayed. The signal is directed 
5 over line 1 8 to the first display monitor 1 2 on which the colored image is displayed. The digital RGB signal is also sent 
over line 20 to the second display monitor 14, which could be located at a location different than the location of the 
first display 12. 

[0031] Prior to applying the RGB signal to the second monitor 1 4, there is interposed an adaptor 22 which receives 
the digital RGB signal and after processing the signal, outputs a signal comprising a new set of digital color values R 1 , 

10 G',B' over line 24 to an input of the second display monitor 14 for each image pixel. 

[0032] The adaptor 22 comprises a receiving and determining element 26 which receives each set of R,G,B values 
and checks them against a plurality of sets of R,G,B values comprising a transform look up table (LUT) 28 which 
correlates R,G,B values to a new set of R'.G'.B' values which when used as an input to the second display monitor 14 
will generate a pixel whose colorimetric value will match that of the pixel generated by the same R,G,B set of color 

15 values in the first display monitor 1 2. 

[0033] If the receiving and determining element 26 identifies a corresponding set of R,G,B values in the LUT 28, it 
outputs the appropriate R'.G'.B 1 values on line 24. If not, the R,G,B input is directed to an interpolator 30 where R'.G 1 , 
B 1 values are derived by interpolation using existing adjacent sets of R,G,B values in the LUT. Preferably, the selection 
of the R,G,B values for the development of the transform LUT 28 is such that values derived by interpolation produce 

20 a visually acceptable match of the displayed pixel in the two monitors 12 and 14. 

[0034] For this system to operate with any degree of success, the development of the transform LUT 28 is critical. 
The values generated by this LUT 28 must indeed lead to values that correctly reproduce colors on the second display 
14 which visually match the colors on the first display 12. Figures 2, 3 and 4 help explain how the LUT 28 is generated. 
[0035] The first step in developing the transform LUT 28 is the creation of two models representing the two displays 

25 whose output is to be matched. Figure 2 shows how this is accomplished. 

[0036] The workstation 10 is used to generate and output over line 16 a set of digital R,G,B color values in a regular 
pattern which is selected to supply color values closely enough spaced so that interpolation of values between existing 
values is reasonably accurate. Preferably, but not necessarily, the output is adjusted so that the color values produce 
displays 32 and 34 on the first and second display monitors 12 and 14, respectively, comprising a plurality of different 

30 color patches on display screens of the display monitors 1 2 and 1 4. As an example, the work station 1 0 may generate 
a plurality of digital RGB values such as shown in the following table I: 



TABLE I 



35 



R: 


0, 


13, 


26, 


51, 


76, 


102, 


128, 


153, 


178, 


204, 


230, 


255. 


G: 


0, 


13, 


26, 


51, 


76, 


102, 


128, 


153, 


178, 


204, 


230, 


255. 


B: 


0, 


13, 


26, 


51, 


76, 


102, 


128, 


153, 


178, 


204, 


230, 


255. 



[0037] Each combination of R,G,B values from this table represents a set (R,G,B) n of color values. In this illustration 
where there are 1728 possible such sets of R,G,B, values, n = 1 to 1728. 

[0038] These 1728 sets of R,G,B, values are supplied to both display monitors 12 and 14 over lines 18 and 20' 
creating a total of 1728 displayed patches for each of the monitors 12 and 14, corresponding to the 1728 different 
combinations of the selected values. All 1728 patches are not necessarily displayed on each monitor screen simulta- 
neously. 

[0039] A colorimetric measuring device 36 is used to read the displayed patches displayed in the first and second 
monitors 12 and 1 4 as shown in Figure 2. The output of the colorimetric device is a colorimetric set of color values for 
each of the patches. In our preferred embodiment, the output of the colorimeter are CI E defined L,a and b sets of color 
values in the Lab color space. Output colorimetric measurements given in other colorimetric color spaces are equally 
acceptable, such as XYZ, Tristimulus CIE defined RGB, etc. Uniform as well as non-uniform color spaces may be used. 
However, values in non-uniform color spaces related to a uniform color space with a known mathematical relation may 
result in increased calculation steps which must be compared with any gained advantage as a result of this selection, 
to determine if the particular choice of color space is justified. (For a description of the different color spaces and related 
terminology, see, in general, DIGITAL IMAGE PROCESSING, 2nd edition, by William Pratt, published by John Wiley 
and Sons, Inc. pages 62-73.) 

[0040] Once the measurement of all patches is completed, two LUTs are compiled, representing the transfer functions 
of the two monitors 12 and 14. The first LUT will consist of the (RGB) n sets of values and the corresponding (Lab) n 
sets of values read off the display on monitor 12 and the second LUT will consist of the same (RGB) n sets of values 
and the corresponding (L'a'b') n sets of values read off the display monitor 14. 
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[0041] Referring to Figure 3, a first and a second model 40 and 42 are built representing the two monitors 12 and 
1 4, respectively. Figure 4 shows a generic model 44 structure, which is used throughout this description in this invention 
to convert color values. The model 44 comprises an LUT 48, which is the particular device derived LUT representing 
the device transfer function. In this example, the LUT 44 for the first model 40 will be the first LUT correlating the (RGB) n 
sets of color values to the (Lab) n sets of values, and the LUT 44 for the second model 42 will be the second LUT 
correlating the same (RGB) n sets of color values to the (L'a'b') n sets of values. 

[0042] In addition to the LUT 48, the model 44 comprises a receiving and determining means 46 similar to, or the 
same as, the receiving and determining means 26 described earlier, and a mathematical interpolating means 50, also 
similar to, or the same as, the previously described interpolating means 30. 

[0043] Figure 3 schematically represents the generation of the transform LUT. Sets of R,G,B color values from a 
preselected plurality of R,G,B sets of color values are used. This plurality of R,G,B sets of values is, preferably, also 
produced in the work station 1 0, and, again preferably, is the same as the sets of (R,G ,B) n values from Table I previously 
used to generate the two LUTs for the two display models 40 and 42. 

[0044] Each (R,G,B) n set of values is directed to the input of the first model 40 over lines 54 and 58. Model 40 
produces an output of (Lab) n values corresponding to the input (R,G,B) n values for this model 40. The same (R,G,B) n 
values are directed to the second monitor model 42 over line 56. Ahead of model 42, there is an adder 60 which 
operates to add to the R n , G n , and B n components of the (R,G,B) n set of values any correction factor dR n , dG n , dB n , 
appearing on line 78. At first, nothing appears on line 78, and the R n , G n and B„ values are applied to the input of 
model 42 unaltered. Model 42 also produces an output set of (L'a'b') n color values corresponding to the input (R,G, 
B) n . This output appears on line 64. 

[0045] The (Lab) n values over line 72 and the (L'a'b') n values over line 64 are next compared in comparator 66 and 
a difference signal [(Lab) n -(L'a'b') n ] is generated and directed over line 74 to correction factor generator 76. Using this 
difference signal, correction factor generator 76 produces correction factors dR^, dG,,^, dB^ (in a manner to be 
explained herein below) which are sent over line 78 and added to the corresponding components of the (R,G,B) n set 
of values in adder 60, resulting in a new set of color values R n (i)'=Fl n +dR n (i), G n(1) '=G n +dG n( i), and B n(1) '=B n +dB n( i) 
on line 62. These new values are again applied to the input of model 42 which produces a new output set of L2'a 2 'b 2 ' 
values on line 64 which is again compared in comparator 66 with the set of (Lab) n values from model 40. A new set 
of correction factors dR n(2) , dG n(2) , and dB n(2) is produced and added to R n(1) ', G n( i)' and B n(1) ' to produce a new set 
of R n ', G n ', and B n ' values. An error "E" defined as: 

E=[(L-L') 2 +(a-a') 2 +(b-b')2] 1/2 

is used to determine when to stop this cycle. Referring to Figure 5, error "E" is compared in comparator 92 with prese- 
lected minimum acceptable limits or a "0" value. If the error is "0", or within the preselected acceptable limits, a signal 
over line 94 switches switch 98 (which is part of adder circuit 60) to feed the R n ',G n ,'and B n ' values which produces 
the minimum over line 80 to the transform LUT compiler 82, rather than to model 42; the next input set of values (R, 
G,B) n+1 is then applied over line 54 to the two models and the whole process is repeated for this new set of values 
and so on. 

[0046] The set of R n 'G n 'B n ' values which produced this zero or acceptable minimum error are sent over line 80 to 
the LUT compiler 82 where the (R,G,B) n and corresponding (R'G'B') n sets of values are used to compile the transform 
LUT 28 so that for each set of (R,G,B) n values there is a corresponding set of (R',G',B',) n values. Means to store the 
transform LUT 28 is included in the compiler 82. This LUT 28 is the the same transform LUT shown in Figure 1 discussed 
earlier. 

[0047] The correction factor generator 76 produces the dR, dG and dB correction factors as follows. The input to the 
correction factor generator 76 is the difference signal produced by the comparator 66, specifically [(Lab) n -(L'a'b') n ]. 
The individual components of this signal are: dL which equals (L-L 1 ), da which equals (a-a') and db which equals (b- 
b'). The correction factors dR, dG and dB, and the difference signals dL, da, and db are related by the following rela- 
tionships: 

(1) (3b/9G)da-(3a/dG)db 
= [ (3b/aG)*(3a/aR) - (3a/3G)*(3b/3R)]dR 



[(3b/3G)*(3a/3B) - (daJdG)*(db/dB)]dB 
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(2) (3a/dG)dL-(aL/3G)da 
= [(da/dG)*(d\JdR) - (3L7dG)*(aa/9R)]dR 
+ [(da/dG)*(dUdB) - (3U3G)*(aa/aB)]dB 



(3) dL = (3L79G)dG + (3UdR)dR + (3L79B)dB. 

10 

[0048] The actual numeric values for 3L/3R, 3L/3G, 3L/3B, etc. for use in solving the above system of equations, are 
obtained by the process schematically illustrated in Figure 5. The input R n ,G n ,B n values are changed by holding two 
the same and incrementing the third by one unit in incrementing module 96. This module comprises a memory for 
temporarily holding the input values of R n , G n , and B n and an arithmetic means for incrementing each of those values 

is by 1 and for applying various combinations of R n , G n , B n , R n +1, G n +1, and B n +1 sets of values to the input of model 42. 
[0049] An example which uses illustrative RGB values, will be used to explain the derivation of the numerical values 
for the partial derivatives needed to solve the above equations, (1 ) (2) and (3). Let the input FL,, G n , and B n values for 
a set of (R,G,B) n values be 100, 100, 100. Assume that when these values are applied to both models 40 and 42, two 
sets of Lab values, (Lab) n and (L'a'b') n are produced such that the error E as previously defined is different than "0" or 

20 a preselected acceptable minimum. Incrementing module 96 produces over line 65 three new color value sets, Rn.Gn, 
(Bn+1), having the numerical values (100, 100, 101) in this illustrative example, Rn,(Gn+1),Bn, having the numerical 
values (100,101,100) and (Rn+1),Gn,Bn, having the numerical values (101,100,100). These new values are applied 
to model 42 over line 63'. The resulting variations in the components of the (L'a'b'),, values produced by model 42 are 
directed to an arithmetic calculator 82 over line 75 and to buffer memory 84. The calculator 82 performs the operations: 

25 (L-L')/R-(R+1), (L-L')/G-(G+1), (L-L')/B-(B+1 ), (a-a')/R-(R+1), (a-a')/G-(G+1), (a-a')/B-(B+1), (b-b')/R-(R+1), b-b'/G- 
(G+1) and b-b'/B-(B+1), to derive the partial derivative values (3L/8R, 8L/3G, etc.) used in equations (1), (2) and (3) 
above. 

[0050] The results of the calculator 82 are supplied over line 83 to another calculating unit 86. The dL, da, and db 
values from comparator 66 are also directed to this calculating unit 86. Unit 86 solves the system of equations (1 ), (2) 

30 and (3) for dR, dB and dG to derive the correction factors. 

[0051] A problem is sometimes encountered in determining a set of values that will produce a best color match on 
two different displays for the same input set of color values. The problem is that a minimum error may be arrived at, 
which may not be an absolute minimum, but only what may be referred to as a local minimum, and there may be other 
combinations of values (R', G', B') n which could produce a smaller even a zero error and a better color match. To test 

35 for this possibility, an optimization operation is performed. A buffer memory 90 holds the input (R,G,B) n values, the 
corresponding (R'.G'.B'Jn values and the error "E n " resulting from the use of this combination. The next set of (R,G, 
B) n+1 values is processed and it, the corresponding (R',G',B,) n+1 values, and resulting error "E n ", are also stored. The 
following (R,G,B) n+2 , (R',G , ,B') n+2 and "E 3 " are stored, etc., until a set of (R,G,B) n+x and (R'G'B') n+x values for which 
the error "E x "=0 is encountered. Once this true "0" error has been identified, a backtracking operation is performed to 

40 determine if a "0" error can be obtained by correcting the set of (R,G,B) n=x .i values using as an input to the model 42 
the (R,G,B) n=x set of values (which is the set for which a "0" error was obtained) while model 40 receives the (R,G, 
B )n-x-i set - Tnis operation will produce a new set (R1 ',G1 ', B1 ') n=x .., of values and a new error "E1 " n=x -i for the in P ut 
values (R.G.B),,^. If "E1" n= x-1 = "0", the (RI'.GI'.BV),,^ set of values is substituted for the (R'G'B 1 ),,^ set of values, 
and the transform LUT28 is recompiled using these (RI'.GV.BI'^x-i values as corresponding to the (R,G,B) n=x _i set. 

45 [0052] If no "0" error is obtained with the above process, then the earlier arrived at result minimum is compared and 
if lower than the newly found minimum, the result stays. Otherwise the R'.G'.B' values giving the least error are used 
in compiling the transform LUT 28 as they will result in a more accurate match. The optimization process is then 
repeated for the next (moving backwards) set of (R,G,B) n=x . 2 input values and sequentially for all (R,G,B) n sets of 
values for n<x always using the values from the set with "E y " = 0, whose n y number is the closest to the tested set n 

so number. When all stored sets of (R,G,B) n values have been tested in this manner, the memory contents of the buffer 
memory 90 are erased and a new set of input (R,G,B) n beginning with the next n, i.e., n=x+1 is applied over line 54 
and the process of compiling the transform LUT 28 continues. 

[0053] For illustration purposes, lets assume that an input set of values (R,G,B) 20 = (100, 100, 100) is applied to 
both the first model 40 and the second model 42 and that this results in a corresponding set of values (R'.G'.B'^q = 
55 (100, 102, 103,) which gives a minimum acceptable error but different from zero. Let us further assume that the next 
set of input values is (R,G,B) 21 and that (R,G,B) 21 = (110, 100, 100). Also let us assume that the corresponding set of 
values (R'.G'.B')^ are (105, 100, 101), and that the error in this instance is zero. 

[0054] When the zero error is detected, the optimization process begins, and no new input sets of (RGB) n values 
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are applied to line 54. 

[0055] The set preceding the set which had the zero error, i.e., set of input values (R,G,B) 20 = (100, 100, 100) is 
applied to the input of model 40. However, the set (R',G',B') 21 = (105, 100, 101) is applied to the input of the second 
model 42 (rather than the set (R,G,B) 20 ) and a new set of corresponding (R1'G1'B1') 20 values is arrived at, i.e., 101, 
5 102, 103, which will either give a zero error or another minimum error. If the result is a zero error or if the new error is 
smaller than the previous, the new (R1'G1 , B1 , ) 2 o values are substituted for the original values (R , G'B , ) 2 o that were 
used earlier in the compilation of transform LUT 28 and will produce a closest match. 

[0056] The process described is best performed in a computer, or if in hardware, under the control of a computer 
suitably programmed to perform the above computations. The diagram of Figure 5 is provided for illustration, rather 
10 than as an actual working embodiment of hardware sufficient to perform the above process, so as to enable one to 
properly program a computer for this task. 

[0057] Referring again to Figure 3, there is shown a block 70 between the output of model 40 and the input to 
comparator 66. This represents an optional color gamut correcting device, which may be used to match the output 
color gamuts of the two display devices. Gamut matching is often often used because the color range that is possible 

is to reproduce with a particular display is usually not identical with the color range possible to reproduce with another. 
Thus, while both displays have a large number of colors that are reproducible by both, there may be certain combina- 
tions of sets of color values that are outside the reproduction capability of one or the other display. While the use of 
the transform LUT 28 in itself inherently provides color gamut matching between the two devices considered, additional 
color gamut matching may be used in the development of the LUT 28, by inserting a gamut matching device 70 as 

20 shown. This device will perform a preselected mathematical manipulation to the (Lab) n set of values so that there will 
not be any sets for which there are no corresponding (L'a'b') n values from the model 42. Such mathematical manipu- 
lation may be as simple as simple truncation or as complicated as data compression. 

[0058] These types of mathematical operations are well known in the art of color matching. For a brief discussion of 
gamut matching and the necessity of color compression, also see "COLORIMETRY FOR ELECTRONIC IMAGING 
25 SYSTEMS, An intensive short course for Scientists and Engineers", Copyright 1 991 , Munsell color science laboratory, 
Rochester Institute of Technology, by R. S. Berns and R. S. Hunter, Chapters 27 and 28. 

[0059] Figures 6, 7 and 8 illustrate the case where the images that are to be matched are images produced in color 
proofers or printers which use a typical four color reproduction system with four color inks, such as Yellow (Y), Magenta 
(M), Cyan (C) and Black (K). 

30 [0060] Referring to Figure 6, as in the previous example, a work station 1 0 generates digital color values representing 
pixels of a color image. Since a four color system is used and the reproducing equipment accepts YMCK signals, the 
work station 10 output is a set of YMCK signals over output line 116. The signals are directed to a first printer 112 over 
line 118 and to a second printer 114 over line 120. Ahead of the input to printer 114 is placed an adaptor 122 which 
again comprises a transform LUT 128, a receiving and determining element 126 and an interpolator 1 30. The LUT 128 

35 is a transform LUT which correlates a plurality of YMCK values to a second plurality of Y'M'C'K' values such that colored 
pixels produced using a set of YMCK values in printer 112 and a corresponding set of Y'M'C'K 1 values from the LUT 
1 28 in printer 1 1 4 will be identical or virtually identical to the standard CIE average observer when viewed under similar 
illumination and surroundings. 

[0061] As for the previous example the interpolator 130 is used to supply corresponding values for YMCK input 
40 values not listed in the transform LUT 128. The transform LUT 128 has been compiled using substantially the same 
process as for the transform LUT 28 used in the previous example. As illustrated in Figure 7, the work station 10 may 
generate a plurality of digital YMCK values such as shown in the following Table II: 



TABLE II 



Y: 


0, 


13, 


26, 


51, 


76, 


102, 


128, 


153, 


178, 


204, 


230, 


255. 


M: 


0, 


13, 


26, 


51, 


76, 


102, 


128, 


153, 


178, 


204, 


230, 


255. 


C: 


0, 


13, 


26, 


51, 


76, 


102, 


128, 


153, 


178, 


204, 


230, 


255. 


K: 


0, 


13, 


26, 


51, 


89, 


128, 


178, 


255. 











SO 

[0062] Each combination of Y,M,C,K, values from this Table II represents a set (Y,M,C,K) n of color values. In this 
illustration, there are 13824 possible such sets of Y,M,C,K, and n = 1 to 13824 values. 

[0063] The 1 3,824 patches resulting from the different color value combinations above are printed by both printers 
1 1 2 and 1 1 4 to provide two sets 1 32 and 1 34 of 1 3,824 patches each. Again as before, a colorimetric measuring device 
5 5 36 is used to read all of the above patches 1 32 and 1 34 and the colorimeter output is used to compile two LUTs 
representing the transfer functions of the two printers 112 and 114. Again it is preferred that the colorimeter output is 
in the Lab color space and the LUTs correlate (Y,M,C,K) n inputs to (L,a,b) n outputs for each of the printers 112 and 
114. These LUTs are used to construct models 140 and 142 using the generic model diagram shown in Figure 4 for 
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the two printers 1 1 2 and 1 1 4, respectively. 

[0064] Figure 8 illustrates the process followed for the compilation of the transform LUT 1 28. This process is similar 
to the process described earlierforthe compilation of the transform LUT 28. The LUT, however, correlates two pluralities 
of color sets, each set having four rather than three variables. Of the four C n ,Y n ,M n ,K n variables, the three C n ,Y n ,M n 
5 variables are independent variables, but the fourth, K n , is not. K n is dependent on the other three from which it is 
derived. The relationship between CYM and K is well known in the graphic arts. 

[0065] Preselected sets of (C,Y,M,K) n values from the work station 10 or elsewhere, are directed to the inputs of 
both model 140 and 142 over lines 154, 162, 156 and 158. It is prefered that the preselected values be the same as 
the values used to compile the LUTs used in the models 140 and 142 and for purposes of this description the same 

10 value combinations resulting using the values in Table II. 

[0066] Each input set of values (C,Y,M,K) n is applied directly to the input of model 140 over lines 154 and 158. The 
same set of values is also applied over line 156 to one input of adder 160. The adder output is directed to switching 
circuit 198. The output of switching circuit 198 is controlled by a signal over line 194 and may be directed to the input 
of model 142 or to transform LUT compiler 182 over lines 162 and 180, respectively. The input set of values is also 

is applied to compiler 1 82 over line 1 55. 

[0067] The adder 1 60 has a second input on which are applied correction factors dC n , dY n and dM n produced in the 
correction factor generator 176. When the (Y,M,C,K) n set of values first appears at the adder 160 input, there are no 
correction factors applied to the adder and the values Y n ,M n ,C n ,K n are applied to the input of model 142 unaltered. 
Model 140 produces an output set of values (L,a,b) n and model 142 produces another set (L',a',b') n . These two sets 

£o are compared in comparator 1 66. For this description, it is assumed that an optional gamut adjustment module 1 70 is 
not used to modify the (L,a,b) n output of model 140 prior to applying it to the comparator 166 over line 172 
[0068] The comparator 166 produces a difference signal [(Lab) n -(L'a'b , ) n ]. This output is directed over line 174 to 
correction factor generator 176. Using this difference signal, correction factor generator 176 produces three correction 
factors dY n(1) , dM n(1) , dC n(1) (in a manner similar to the manner in which correction factor generator 76 produced the 

25 dR, dG, and dB factors as previously explained). The correction factors are sent over line 178 and added to the cor- 
responding components of the (Y,M,C,K) n set of values in adder 60, resulting in a new set of color values 

Y 'n(1)= Y n +dY n(1). M 'n(1)= M n +dM n(1). and C 'n(1)= C n +dC n(1) on line 162 

[0069] The fourth value, K n is stored in element 1 77 which comprises a buffer memory and a means to modify the 
stored value K n by predetermined increments. At this time K n remains constant and these new Y' n(1) , M' n(1) , C' n(1 } , and 
30 K n values are again applied to the input of model 1 42 which produces a new output set of Lg'ag'^' values on line 1 64 
which is again compared in comparator 166 with the set of (Lab) n values from model 140. A new set of correction 
factors dY n(2) , dM n(2) , and dC n(2) is produced and added to Y' n(1) , M' n(1) and C' n(1) to produce a new set of Y n ', M n ', 
and C n ' values. Again, an error "E" is defined as: 

35 9 2 1/2 

E=[(L-L') 2 +(a-a') 2 +(b-b')2] 1/2 

[0070] The correction factor generator 176 also includes means to calculate, compare and store this error E n and 
the corresponding (YMCK)n and (Y'M'C'K)n set of values that resulted in this error E n . 

40 [0071] The above process is repeated continuously correcting the original input values Y n ,M n ,C n , with K n held con- 
stant until the error E is zero or a predetermined acceptable minimum. If a zero error is found, then the two sets of 
values produce a perfect match and the search terminates. Switch 1 98 is actuated and the two sets of values are used 
in the compiler 182 to compile the transform LUT 128. If only a minimum is arrived at it is again stored in the buffer 
memory and K n is incremented in unit 177 and the process repeated in search of another set of (Y"M"C"K') n values 

45 yielding another minimum. This new minimum is compared with the previously found one. The (Y'M'C'K),, values or 
the values (Y"M"C"K') n which produced the best minimum are kept. 

[0072] The correction factor generator in this instance solves a system of differential equations that is of the same 
format as previously discussed. The equations only differ in terms of variable designation and are as follows: 

50 (1a) (ab/3Y)da-(aa/3Y)db 

= [(3b/aY)*(3a/dM) - (da/3Y)*(ab/3M)]dM 
+ [(3b/aY)*(3a/3C) - (9a/3Y)*(db/3C)]dC 



(2a) (aa/3Y)dL-(3L/aY)da 
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= [(3a/3Y)*(3L/3M) - (3L/3Y)*(3a/dM)]dM 
+ [(9a/3Y)*(dL/aC) - (3L/dC)*(3a/3C)]dC 



(3a) dL = (3L/dY)dY + (3L/3M)dM + (3L/3C)dG. 

[0073] The partial derivatives are calculated in the same manner as the partial derivatives were calculated for the 

10 R,G,B[ case, by incrementing the Y n ,M n ,C n values by one, as explained earlier with the help of Figure 5. 

[0074] The issue of whether a minimum error E n is a true minimum error or a local minimum, meaning that there 
may be another smaller error resulting from a different combination of (C',Y',M',K') n values which could be found by 
starting with a different input set to the second model 142, rather than the (Y,M,C,K) n set that was used is again a 
problem. The same testing and optimization process described earlier can be used in this case too, to test and discover 

15 whether there may be other better combinations of (Y',M',C',K') n values, by holding in a buffer memory all the consec- 
utive sets of (Y,M,C,K) n , the corresponding (Y'M'C'K) n sets and error E n as described earlier. However, it has been 
found that more accurate results may be obtained by following a somewhat more calculation intensive process. This 
process utilizes two input sets of values which have resulted in a zero error as inputs to obtain two optimized sets of 
(Y1',M1',C1',K1') n and (Y2 , ,M2',C2',K2 , ) n values in addition to the (Y'.M'.C'.K'Jn set and then selecting the set that 

20 produces the best match, i.e., least error, from the three sets. 

[0075] Because of the presence of the fourth color K, which is not an independent color as the Y,M,C, and which by 
its dependency on the other three can introduce a number of color combinations resulting in matching colors, the 
existence of two different sets of values both producing a perfect match is also possible. Further, the optimization 
process must be able to rationally select the set that will produce the better overall match if used in the compilation of 

25 the transform LUT 1 28. When two sets produce perfect matches, then the one that produces better interpolation results 
for input values that are not present in the LUT will be selected. 

[0076] The optimization process adopted for this case, is based on selecting a predetermined fixed number "x" of 
sets of values (Y,M,C,K) n , the corresponding (Y',C',M',K') n and the associated error E nx and, then, proceeding with the 
optimization calculation using zero error resulting set values from either side of a target set to derive two sets of possibly 

30 lesser error resulting values in addition to the set originally found. Thus, the input set (Y,M,C,K) n8 , for instance, for 
which the set of (Y'C'lv1'K') n8 gave us originally a minimum but non zero error, will be optimized using preceding and 
following values, to obtain (Y1\ C1\ M1',K1') n8 values with data from a set (Y',C',M',K') ny and (Y'2,M'2,C'2,K'2) n8 using 
data from a set (Y , ,C',M',K') n2 where y>8 and z<8 and where both the "ny" and "nz" sets are the sets which produced 
a zero error closest to the set n=8. 

35 [0077] The following illustrative example is provided to explain the optimization process based on a fixed number of 
value sets discussed above. 

[0078] Assume that the transform table 1 28 is to be constructed using the (YMCK) n sets of color values from Table 
II. Let us further assume that we are at a stage in the process where the sets of (YMCK) n shown in Table III have been 
sequentially introduced over line 1 54 to model 1 40, 1 ine 1 56 to adder 1 60 and line 1 55 to LUT compiler 1 28. As a result, 
40 corresponding sets of (Y'M'C'K') of values and associated error En values have been obtained. The errors En at this 
time have not been tested for optimum error condition, for errors different than zero. The number of sets used for 
optimization will be limited to four (that is x=4), to make the illustration simpler. The four sets of values and associated 
error values held in memory at the beginning of the optimization are shown in Table III. 



Table III 



so 



Set* 


(Y,M,C,K) 


(Y',M',C',K) 


E 


n 


76, 


76, 


76, 


89 


73, 


72, 


78, 


90 


0 


n+1 


76, 


76, 


102, 


89 


73, 


73, 


100, 


89 


0.1 


n+2 


76, 


76, 


128, 


89 


73, 


75, 


131, 


91 


0 


n+3 


76, 


76, 


153, 


89 


73, 


71, 


152, 


93 


0 



[0079] For input sets n, n+2 and n+3 the error is zero. Thus, we have perfect match and can accept the (Y,M,C,K) 
and corresponding (Y',M',C',K') values in the compilation of the transfer table 128. However, we need to check the 
results of set n+1 to see if a better match may be available. 

[0080] We proceed as we did earlier in the optimization process discussed for the (R,G,B) example. However, the 
process is done with two sets of input values to model 142, the sets that gave a zero match located on either side of 
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the set that is optimized. In this example, the set that is optimized is set n+1 and the two sets, that gave zero error on 
the two sides of set n+1 , are the n set and the n+2 set. We apply the (Y,M,C,K) n+1 (i.e., 76, 76, 102, 89) set of values 
again to model 140, but we apply to model 142 through the adder 160 first the (Y',M',C',K') n (i.e., 73, 72, 78, 90) set 
of values which is the set corresponding to the "n" input, the set which gave a zero error in matching the input values 
5 for "n". 

[0081] The (Lab) n+1 and (L'a'b') n outputs from the two models are compared in comparator 166 and new correction 
factors dY' n , dM' n , dC' n , are generated which result in a new set, Y'1 n+1 , M'1 n+1 , C'1 n+1 , K'1 n+1 of values when added 
to the (Y\ M',C, K') n in adder 1 60. As before, a new set of L'a'b' values is produced and again compared in the comparator 
166, and the process repeated, including the K incrementing step, until a minimum error value is found again. Assume 
10 these new values to be Y'1 n+ i=73, M'1 n+ i=74, C1 n+1 =1 00 and K'1 n+1 =88. Let us further assume that this set of values 
results in a perfect match and a zero error. We temporarily store these results. 

[0082] Next we repeat the above process using this time the set of (Y'M'C'K') n+2 values as an input to the adder 1 60, 
while the (Y,M,C,K) n+1 set continues to be applied to model 140. This time we find that another set of 
(Y^.M^.C^.K^,),,^ values, i.e., Y'2 n+1 =71, M'2 n+1 =76, C'2 n+1 =100 and K'2 n+1 =91, where the error is again zero. 
is These values are also stored. These result are shown in Table IV: 



TABLE IV 



Set# 


(Y,M,C,K) 




(Y'.M'.C'.K') 


E 




76, 


76, 


76, 


89 




73, 


72, 


78, 


90 


0 


n+1 


76, 


76, 


102, 


89 


(a) 


73, 


73, 


100, 


89 


0.1 












(b) 


73, 


74, 


100, 


88 


0 












(c) 


71, 


76, 


100, 


91 


0 


n+2 


76, 


76, 


128, 


89 




73, 


75, 


131, 


91 


0 


n+3 


76, 


76, 


153, 


89 




73, 


71, 


152, 


93 


0 



[0083] The question now becomes which set of values (Y',M',C',K') n+1 should we keep in compiling ourtable 128 for 
the input set (Y,M,C,K) n+1 . Obviously set (a) is the worst match, and we can discard it. Of the other two, it is preferred 
to keep the set that will provide the most accurate interpolation results, when interpolation involving those values must 
be used. One way to do this is to select the set that has values best distributed between the preceding "n" set and the 
following "n+2" set. A way to make this selection is as follows. Let Table V represent the two possibilities of selection. 



TABLE V 



40 



Set# 




Possibility #1 






Possibility #2 






Input 






(Y'.M'.C'.K) 






(Y'.M'.C'.K') 






(Y,M,C,K) 


n 


73, 


72, 


78, 


90 


73, 


72, 


78, 


90 


76, 


76, 


76, 


89 


n+1(b) 


73, 


74, 


100, 


90 










76, 


76, 


102, 


89 


n+1(c) 










71, 


76, 


100, 


88 


76, 


76, 


102, 


89 


n+2 


73, 


75, 


131, 


91 


73, 


75, 


131, 


91 


76, 


76, 


128, 


89 



[0084] We next compare the point spread deviation between input (Y,M,C,K) sets and the corresponding (Y'.M'.C, 
K) sets. For instance, comparing the inputs between sets n and n+1, we have: Y n+1 -Y n =0, M n+1 -M n =0, C n+1 -C n =26, 

45 K n+1 -K n =0 and comparing the inputs between sets n+1 and n+2, we have: Y n+2 -Y n+1 =:0, M n+2 -M n+1 =0, C n+2 -C n+1 =26, 
K n+2 -K n+1 =0. Next, we do the same using the values for Y'.M'.C'.K' from the first selection column (labelled "Possibility 
#1), and we get: Y' n+1 -Y' n =0, M' n+1 -M' n =2, C' n+1 -C' n =22, K' n+1 -K' n =-1 . Thus, we see that for a zero change in the Y 
input values, there is a zero change in the Y" values. For zero change in the M input values, we have a change of +2 
points in the M' value. For a 26 point change in the C value input, the C value changes by 22 points, or a short fall of 

so -4. And for a zero input change in the K value we have a -1 point change in the K'. 

[0085] Similarly the differences observed between the n+2 and n+1 sets, yields a difference for Y=0 points, M=+1,, 
C=+5 and K=+2. 

[0086] We will refer to these differences in the values as the deviation from the ideal situation where the increments 
of change between input values are the same as the increments of change between corresponding values, that is 
55 vibere{C n -C n _- l )=(C' n -C' n .i) . The average of the deviation magnitude AD n for the deviation between sets n and n+1 for 
Y'.M'.C'.K' from Table V selection #1 column is AD 1 =7/4 or 1 .75; for sets n+1 and n+2, again for selection #1 the average 
is AD 2 =2. 

[0087] For values from column selection #2 in Table V, the same process yields AD 3 =2.75 for sets n and n+1 , and 



EP 0 604 755 B1 



2.25 for AD 4 for sets n+1 and n+2. Thus, after comparing the average deviations for the two selections, set n+1(b) in 
selection #1 is seen to yield set values that have magnitudes which, when taken as adjacent pairs of values with the 
preceding or following sets of values define color spans that on the average deviate less from the spans defined by 
the corresponding input color pairs than the values obtained from the sets of values from selection #2. Therefore, the 
5 set of (Y',M',C',K') n+1(b) values is used in the transform table compilation. If set (a) could not be discarded as obviously 
the worst match, it could be compared to sets (b) and (c) like set (b) was compared to set (c) above to determine which 
set is the best match. 

[0088] Figure 9 shows a situation where image color matching is desired in images produced with a color printer 214 
which accepts RGB values (through a converter 235 that converts RGB inputs to YMCA using a mathematical rela- 

10 tionship that the printer manufacturer has determined as an accurate color conversion) and a CRT monitor 212. For 
this arrangement, as shown in Figure 9, the output of workstation 10 is an RGB signal over line 216. It is used to drive 
the CRT display 212 directly and the color printer 214 through a 3 to 4 color conversion unit 225 which converts RGB 
inputs to YMCK outputs; an adaptor 222 comprising a transform LUT 228 constructed in accordance with this invention 
is placed ahead of conversion unit 225. 

rs [0089] This LUT 228 is compiled using a process similar to the process illustrated in Figure 8. However, in generating 
the model 44 for the printer 21 4, the converter unit 225 is included with the color printer 21 4 as an integral part thereof. 
Th us, this model 44 which is schematically shown in generic format in Figure 1 1 has an input 45' which leads to converter 
47. The output of converter 47 over line 45 is directed to the input of a receiving and determining element 46 and 
therefrom in a similar arrangement to that shown in the generic model of Figure 4 produces an output 52 either directly 

20 through LUT 48, or by interpolation through interpolator 50. 

[0090] In the arrangement shown in Figure 1 0, the output of the work station 1 0 is a YMCK signal which drives directly 
color proofer 312 over lines 316 and 318. The YMCK output is also used to display an image on a CRT monitor 314. 
This monitor 314 is able to accept YMCK input values because of a 4 to 3 color conversion that occurs in a color 
converter 325 prior to being used by the monitor 314. To obtain more accurate image color matching, an adaptor 322 

25 which includes a receiving and determining element 326 and a transform LUT 328 for converting YMCK input values 
to Y'M'C'K' is placed ahead of the converter 325. As before an interpolator 330 is included in the adaptor to derive 
values not present in the LUT 328. The Y'M'C'K' values are directed to a CRT display 314 over line 324 through existing 
converter 325. 

[0091] The transform LUT 328 is also generated using the process described earlier, wherein two models 44 are 

30 created representing the two display devices 312 and 314. The CRT display device model 314, however, as in the 
previous example, includes the 4 to 3 converter 325 and is also similar to the generic model illustrated in Figure 11. 
The LUT 328 is again compiled using the process illustrated in Figure 8. The transform LUT 328 converts input YMCK 
values to output Y'M'C'K' which when sent to converter 325 generate RGB values such that when displayed on the 
CRT display 314 will produce an image having the same or substantially the same color appearance as the image 

35 produced by the color proofer 31 2 using the original YMCK values directly. 

[0092] In discussing this invention there has been constant reference to "interpolation". Interpolation leading to de- 
rived output sets of color image values based on input sets of image values is used in performing the process of this 
invention. It is evident that if interpolation is not used, the LUTs employed both in the models and in the adaptors must 
supply values for all possible combinations of input color values. In a digital system operating in an 8bit domain, for a 

40 three variable system, i.e., RGB, the possible input value combinations are 256*256*256 or 1 6,777,21 6 combinations. 
To develop the necessary LUTs, each of these combinations would have to be printed as a patch, in accordance with 
this process, measured, compared, and optimized, clearly an impractical process because of the shear magnitude of 
numbers and calculations involved. Thus, limited numbers of sets are selected as illustrated in Tables I and II and 
interpolation is used for intermediate sets of input color values. 

45 [0093] The particular interpolation process used may vary. Depending on the accuracy one wishes in the system 
more or less calculation intensive interpolations may be used. In practicing this invention, it is preferred that a less 
calculation intensive interpolation method is used in the model interpolators, while a more sophisticated method is 
used in the adaptor interpolators. Examples of both types are given below. However, any of the interpolation methods 
described herein below, as well as other interpolation formulae, such as quadratic interpolation, can be used for either 

so the adaptor interpolators, or the model interpolators, the final choice being based on desired accuracy of the matching 
process, and the availability of resources. This invention is not restricted to the use of one interpolator for adaptors 22, 
1 22, 222, 322 and a different one for the models 40, 42, 1 40 and 1 42 even though such arrangement is both used and 
preferred. 

[0094] The preferred interpolation performed in interpolators 30, 1 30, 230 and 330 shown in the adaptors shown in 
55 Figures 1,6, 9 and 10, respectively, that is the interpolator 30, 130, 230 or 330 that accompanies a transform LUT 28, 
1 28, 228 or 328, respectively, is tetrahedral interpolation performed through a programmed computer. A pair of appro- 
priate programs for programming such a computer to perform this interpolation for three and four variables are included 
with this application as Appendices "A" and "B" which appear immediately before the Claims. 
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[0095] Linear interpolation, tri-linear or quadralinear depending on the particular color space in which the color values 
exist, is the preferred choice of interpolation in models 40, 42, 140 and 142 shown in Figures 3 and 8, used in deriving 
the transform LUTs 28, 128, 228, 328, as it is simpler, and often adequate to produce good results. The algorithm for 
this type of interpolation as applied in this invention is next explained for the more complex case of a four to three 
5 transform, i.e., from (YMCK) input values to (Lab) as would be used in models 140 or 142. The three to three (RGB 
color space) and four to four (YMCK color space) transform cases are similar. 

[0096] The numerical values for Y,M,C,K given in Table II represent data points along a system of axes Y, M, C, and 
K. In the case of three variables, i.e., R,G,B color values, this may be an orthogonal axis system that may be easily 
visualized. With four variables visualization is very difficult, however, such actual visualization is not necessary to 

70 perform the interpolation. 

[0097] The data points representing predetermined color values of a variable along the axis for this variable are 
separated by a specified number of units representing a distance measured along this axis. Thus, the data points 
Y=102 and Y=128 are separated by a distance of 26 units along the Y axis, while the data points M-128 and M-153 
are separated by 25 units along the M axis. An input value may fall either on a data point, or between data points. 

is When it falls on a data point the LUT provides a corresponding value. If the input value fall between data points, then 
an output is calculated as follows. 

[0098] Let F| be the ratio of the distance of the input value from a preceding data point along an axis to the distance 
between that data point and the next consecutive data point along this axis. Thus, if the color value set has four 
components, (i=0, 1 , 2, or 3, for the Y, M, C, K components) and Y n is a data point along the Y axis, Y n+1 is the next 
20 data point along this axis, and Y-, is the Y component value of a color input set of values whose numerical value lies 
between Y n and Y n+1 , then: F F (Y1-Y n )/(Y n+1 -Y n ) for kO. 

[0099] We also define a quantity Qi=1-Fi, and the coefficients tj (where j=1, ... 16) as follows: 



t 1 =F 0 *F 1 *F 2 *F 3 


t 6 =Q 0 *CVF 2 *F3 


I^Fq'F^CVQs 


t 2 =Q 0 *F 1 *F 2 *F 3 




t 12 =F 0 *QrQ 2 *Q 3 


t 3 =F 0 *Q/F 2 *F 3 




t 13 =Q 0 *F 1 *Q 2 *Q 3 


t 4 =F 0 *F 1 *Q 2 *F 3 


tsFrVQi*Q 2 *F 3 


^ 4=00*0! *F 2 *Q 3 


t 5 =F 0 *F 1 *F 2 *Q 3 


tiO=Fo*Qi* F 2*Q3 


t 15 =Q 0 *Q 1 *Q 2 *F 3 


ti6=Qo*QrQ 2 *Q 3 







[0100] Let L( Yn Mn cn Kn). a(Yn,Mn,Cn,Kn)> b (Yn,Mn,Cn,Kn)> represent the Lab values from the LUT for the model 140 for 
instance, for an input Y n ,M n ,C n ,k n ,' L( Yn+1iMniCn , Kn ). a( Y n+i,Mn,c n> Kn)> b( Yn+ i,Mn,cn,Kn) the Lab values for an input Y n+1 , 

M n , C n , and K n , L( Yn , Mn+1 , C n,Kn)> a( Yn , M n+1,Cn,Kn). b( Yn , M n + 1,Cn,Kn)> the ValU<3S fr ° m the LUT f ° f a " ' n P Ut Y "' M n+T C n 

3S and K n , and so on for all 16 combinations of input values (Y, M, C, K) n and (Y, M, C, K) n+1 . The Lab values L v a v and 
b., for an input set of Y1 ,M1 ,C1 ,K1 color values falling between the n and n+1 sets of values in the LUT, are given by 
the following equations: 

L 1 = l 1 * L (Yn+1,Mn+1,Cn+1,Kn+1) + t2 * L (Yn,Mn+1,Cn+1,Kn+1) 

40 

+ t3 * L ( Yn +l ,Mn,Cn+1 ,Kn+1 ) + 14 * L (yii+1 ,Mn+1 ,Cn,Kn+1 ) 
+ t5 * L( Yn+ i ; Mn+1,Cn+1,Kn) + t6 * L (Yn,Mn,Cn+1 ,Kn+1) 

45 + 17 * L( Yn Mn+1 Cn Kn+1 , + 18 * L( Yn Mn+1 Cn+1 Kn) 

+ t9 * L ( Yn+1 , M n,Cn,Kn + 1) +t1 ° * L <Yn + 1 ,Mn,Cn + 1 ,Kn) 
+ t11 * L ( Yn+1, Mn+1 ,Cn,Kn) +t12 * L (Yn+1,Mn,Cn,Kn) 

SO 

+ t1 3 * L (Yn,Mn+1 ,Cn,Kn) + 11 4 * L (Yn,Mn,Cn+1 ,Kn) 
+ t15 * L( Yn>Mn C n,Kn+1) + 116 * L (Yn,Mn,Cn,Kn)- 



55 

Similarly, 
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a 1 =t, *a( Y n + 1,Mn+1,Cn+1,Kn+1) + 12 * a (Yn,Mn+1, Cn+1, Kn+1) 
+ t3 * a (Yn+1 ,Mn,Cn+1 ,Kn+1 ) + t4 * a (yn+1 ,Mn+1 ,Cn,Kn+1) 
+ t5 * a( Y n+1,Mn+1,Cn+1,Kn) + 16 * a (Yn,Mn,Cn+1 ,Kn+1) 
+ t7 * a( Yn ,Mn+1 ,Cn,Kn+1 ) + t8 * a (yn,Mn+1 ,Cn+1 ,Kn) 
10 + t9 * a( Yn+1 iMn C n,Kn+1) + 11 0 * a ( Yn +1 ,Mn,Cn+1 ,Kn) 

+ f 1 * a ( Yn+ i ,Mn+1 ,Cn,Kn) + t1 2 * a ( Y n+1 ,Mn,Cn,Kn) 
+ t1 3 * a ( YnMn+1 ,Cn,Kn) + t1 4 * a ( Yn ,Mn,Cn+1 ,Kn) 
+ t15 * a ( YniMn ,Cn,Kn+1) + 116 * a ( Yn ,Mn,Cn,Kn)- 

And, 

b 1 =t 1 *b( Yn+1 Mn+1 Cn+1 Kn+1) +t2 *b( Yn Mn+1 Cn+1 Kn+1) 
+ 13 * b( Yn+1 Mn Cn+1 Kn+1 j + 14 * b( Yn+1 Mn+1 Cn Kn+1 ) 

25 + 15 * b (Yn+1, Mn+1, Cn+1 ,Kn) + 16 * b (Yn,Mn,Cn+1 ,Kn+1) 

+ *7 * b ( Yn ,Mn+1,Cn,Kn+1) + t8 * b (Yn,Mn+1 ,Cn+1 ,Kn) 
+ t9 * b <Yn+1 ,Mn,Cn,Kn+1) + t1 0 * b (Yn+1 ,Mn,Cn+1 ,Kn) 
+ t11 * b (Yn+1,Mn+1,Cn,Kn) + 112 * b (Yn+1,Mn,Cn,Kn) 
+ t1 3 * b ( Y n,Mn+1 ,Cn,Kn) + 11 4 * b ( Y n,Mn,Cn+1 ,Kn) 
3g + t15 * b ( Y n,Mn,Cn,Kn+1) + 116 * b (Yn,Mn,Cn,Kn)' 

[0101] Those skilled in the art, having the benefit of the teachings of the present invention as herein above set forth, 
can effect numerous modifications thereto. These modifications are to be construed as being encompassed within the 
scope of the present invention as set forth in the appended claims. 

40 



so 
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APPFNTVCX A 



/************************************************************** ****** 

(c) Copyright 1992 
E. I. DuPont do Nemours and Company (inc.) 
Wilmington, Delaware 19898. 

tr3us_t.c 

A function fox tetrahedral type highspeed interpolation 
for performing 3 to 3 or 3 to n color coordinate 
conversions for the case of uniform sampling along 
coordinate axes. The vertices of the cube are 
numbered sequentially in a natural binary order 
(x3 x2 xl) . 

C. B. Chittineni 

June 1, 1992 



tr3us_t(n,ns,xl,x2,x3,l,x,y, z,val) 
no of outputs (int short) 

total number of data sample points per output (int short) 
data point with in the cube along xl axis (int short) 
data point with in the cube along x2 axis (int short) 
data point with in the cube along x3 axis (int short) 
length of cube (int short) 

data vector at current node (least rapidly varying) 
(unsigned char) 

data vector at current node (second least rapidly 
varying) (unsigned char) 

data vector at current node (most rapidly varying) 
(unsigned char) 
val- pointer to interpolated values (n) (int short) 

Modification: (Chittineni: July 9, 1992)Changed the vertex numbering 
as as (n) -> (x y z) (0 - 1 1 1, 1 - 0 1 1, 2 - 1 0 1, 

3-00 1, 4-11 0, 5-01 0, 6-10 0, 
7-00 0) (Liang's data files format). 

A***************,***************************************************** i 

40 tr3us_t(n, ns, xl, x2, x3, 1, x, y, z, new) 

short new[3], n, ns, xl, x2, x3, 1; 
unsigned char *x, *y, *z; 
{ 

int short wO, wl, w2, w3, a, b, c, d; 
45 int itm, j, jl, dt; 

int j2, j3, j4; 

unsigned char *x0, *y0, *z0; 

itm-0; 

if (x2 >- xl) { 
so if (x3 >- x2) { 

itm-3; 
}else if (x3 >- xl) { 
itm-5; 
}else{ 

itm-4; 

ss } 

)else if (x2>-x3) { 



Module : 
Purpose: 



Author : 
Date Create: 

Usage : 

where n ■ 



x2 
x3 
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itm-0; 
}else if (x3 >- xl) { 
itm-2; 
}else{ 

itm-1; 

} 

switch (itm) { 
case 0: 

wO-(l-xl); wl-(xl-x2); v2-(x2-x3); w3-x3; 
n«w[03-(wO*x[7]+wl*xr3]+w2*x[l]+w3*xrO])/l; 
newtl]-(w0*y[7]+wl*y{3]4w2*yll)+w3*yr0])/l; 
new[2)-<w0*z l7]+wl*z {3]+w2*z [l]+w3** [OJ ) /In- 
break; 

case 1: 

wO-(l-xl); wl-(xl-x3); w2-(x3-x2); w3-x2; 
new[0]-(wO*xI7]+wl*x[3J+w2*x[2J+w3*x[0])/l; 
new t 1 ) - <wO*y (7 J +wl*y [3] +w2*y [ 2 ] +w3*y [ 0 ] ) /l ; 
newt2]-(w0*z[7]+wl*zt3]+w2*z[2]+w3*z[0])/l; 
break; 

case 2: 

w0-(l-x3); wl-(x3-xl); w2-(xl-x2); w3-x2; 

new[0]-(w0*x[7]+wl*x[6]+w2*x[2]+w3*x{0])/l; 

new[l]-(wO*yl7)+wl*yl6]-^w2*yt2]+w3*ytO])/l; 

new[2]-(w0*zl7]+wl*z[6J+w2*zt2J+w3*zl0])/l; 

break; 

case 3: 

w0-{l-x3); wl-(x3-x2); w2-(x2-xX); w3-xl; 

newI0J-<w0*x[7]+wl*x[6]+w2*x[4]+w3*x[0])/l; 

new[lj-(w0*y[7]+wl*y[6]+w2*y[4]+w3*y[0])/l; 

newt2]-(w0*z[7]+wl*z{6]+w2*z[4)+w3*z[0])/l; 

break; 

case 4 : 

w0-(l-x2); wl-(x2-xl); w2-(xl-x3); w3-x3; 
newt0]-(w0*x(7}+wl*xl5)+w2*xll]+w3*x[0])/l; 
newtl]-<wO*y{7]+wl*y[5J+w2*y[l]+w3*y(0])/l; 
newfr-2]- (wO*z [7] +wl*z (5J+w2*z [1] +w3*z [0] ) /l; 
break; 

case 5: 

w0-(l-x2); wl-(x2-x3); w2-(x3-xl); w3-xl; 
new(0]-{wO*x[7]+wl*xI5]+w2*x[4]+w3*xtO])/l; 
new[l]-(w0*yf7]+wl*y[5]+w2*yt4]+w3*y[0])/l; 
new[2]-(w0*z [7]+wl*z [5]+w2*z [4]+w3*z [0] ) /l; 
break; 

} 

return; 
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APPENDIX B 



(c) Copyright 1992 
E. I. DuPont de Nemours and Company (Inc.) 
Wilmington, Delaware 19898. 

Module: tr4us_t.c 

Purpose: A function for tetrahedral type highspeed interpolation 

for performing 4 to 4 or 4 to n color coordinate 
conversions for the case of uniform sampling along 
coordinate axes. The vertices of the hypercube are 
numbered sequentially in a binary order (xl x2 x3 x4) . 
It is to be noted that the ordering is important 
in arranging the sampled data and in passing the 
arguments to the function. 

Author: C. B. Chittineni 

Date Create: June 1, 1992 

Usage : tr4us_t (n, ns,xl,x2, x3, x4, 1, cbl, mbl, ybl, kbl, val) 

where n - no of outputs (int short) 

ns - total number of data sample points per output (int) 
xl - data point with in the hypercube along xl axis (int short) 
x2 - data point with in the hypercube along x2 axis (int short) 
x3 - data point with in the hypercube along x3 axis (int short) 
x4 - data point with in the hypercube along x4 axis (int short) 
1 - length of hypercube (int short) 

cbl- data vector at current node (least rapidly varying) 
(unsigned char) 

mbl- data vector at current node (second least rapidly varying) 
(unsigned char) 

ybl- data vector at current node (third least rapidly varying) 

(unsigned char) 
kbl- data vector at current node (most rapidly varying) 

(unsigned char) 
val- pointer to interpolated values (n) (int short) 

Modification: (Chittineni: July 9, 1992) Changed the vertex numbering 
as (n) -> (cbl mbl ybl kbl) (0-1111, 1-0111, 
2-1011, 3-110 1, 4-1110, 5-0011, 
6-010 1, 7-0110, 8-100 1, 9-1010, 
10- 1 1 0 0,11 - 1 0 0 0,12 - 0 1 0 0,13 - 0 0 1 0, 
14- 0 0 0 1,15 - 0 0 0 0) (Liang's data format) . 

******************************************************************* 

♦include <stdio.h> 
♦include <math.h> 

tr4us t(n, ns, xl, x2, x3, x4, 1, cbl, mbl, ybl, kbl, new) 
int snort n, xl, x2, x3, x4, 1/ 
int short *new; 
int ns; 

unsigned char *cbl, *mbl, *ybl, *kbl; 



< 



int short wO, wl, w2, w3, w4, a, b, c, 
int itm, j, jl, j2; 



/* Locate the tetrahedran */ 
if (xl >- x2) { 
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if 1x2 >- x3){ 

if (x3 >- x4){ 

itm-1; 
}else if (xl >- x4){ 

if (x2 >- x4){ 

itm-2; 
}else{ 

itm-3; 

} 

>else { 

itm-19; 

} 

Jelse if (x3 >- x4) { 

if <xl>-x4){ 

if <x2>-x4){ 

if (xl>-x3){ 

itxn-5; 
}else{ 

itra-13; 

}else if (xl >- x3) { 
itm-6; 
)else{ 

itro-14; 

} 

}else{ 

itro-15; 

} 

}else if (xl >- x4) { 

itm-4; 
Jelse if (xl >- x3) { 

itm-20; 
}else { 

itra-21; 

} 

}else if (x2 >- x3) { 

if (x3 >- x4) { 

if(xl >- x4){ 

if (xl >- x3) { 

itm-7; 
}else { 

itm-11; 

} 

}else{ 

itm-12; 

} 

}else if (xl >- x4) { 
itm-8; 
}else if (x2 >- x4){ 

if (xl >- x3) { 

itm-9; 
}else{ 

itra-10; 

} 

}else if (xl >- x3) { 
itm-2 3; 
}else{ 

itm-2 4; 

} 

Jelse if (x3 >- x4) { 

if(xl >- x4)( 

itm-17; 
)else if (x2 >- x4) { 
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ittt-18; 
}else{ 

itro-16; 

) 

)else< 

itm-22; 

} 

fprintf (stderr, "%d \n", itra);*/ 
switch (itm) { 
case 1: 

wO-(l-xl); wl-(xl-x2); w2-(x2-x3) ; w3-(x3-x4); w4-x4; 

new [0] - (wO*ebl [15]*wl*cbl til] +w2*cbl [10 J +w3*cbl [4)+w4*cbl [ 0] ) /l; 

new [1] - (w0*mbl [15J +wl*mbl 111] +w2*nfcl 110] +w3*mbl [4]+w4*mbl 10] ) l\; 

new [2] - (w0*ybl 115] +wl*ybl [11] +w2*ybl [ 10 J +w3*ybl [4] +w4*ybl [ 0] ) /l; 

new [3] - {w0*kbl 115] +wl*kbl [11] +w2*Jebl f 10 ] +w3*kbl [ 4 ] +w4*kbl [ 0 J ) /l ; 

break; 

case 2: 

wO-(l-xl); wl-(xl-x2); w2-(x2-x4); w3-(x4-x3); w4-x3; 

new [0] - (w0*cbl [15] +wl*cbl [11] +w2*cbl [10] +w3*cbl [3] +w4*cbl [0] ) /l; 

new [1 ] - (w0*mbl [15] +wl*mbl [11] +w2*wbl [10] +w3*mbl [3] +w4*mbl [0] ) /It 

new [ 2 ] - ( w0*ybl [ 1 5 ] +wl *ybl [ 11 ] +w2*ybl [ 10 ] +w3*ybl [ 3] +w4 *ybl f 0 ] ) /l ; 

new [3]- (w0*kbl [15] +wl*kbl [11] +w2*kbl [10] +w3*kbl [3J+w4*kbl [0] ) /l; 

break; 

case 3: 

wO-U-xl); wl-(xl-x4); w2-(x4-x2); w3-(x2-x3); w4-x3; 
new(0]-(w0*cbl[15]+wl*cbl[ll]+w2*cbl[8]+w3*C!bl[3]+w4*cbl[0]>/l; 
new [ 1 ] - ( w0*mbl [ 15 ] +wl*nibl [ 11 ] +w2*mbl [ 8 ] +w3*mbl [3] +w4 *mbl [ 0 ] ) /l ; 
new [2] - (w0*ybl [ 15] +wl*ybl [11] +w2*ybl [8 ] +w3*ybl [3] +w4*ybl [0] ) /l; 
new ( 3 ] - ( wO *kbl [ 1 5 ] +wl *kbl [ 1 1 ] +w2 *kbl [ 8 ] +w3 *kbl [ 3 ] +w4 * kbl ( 0 ] > / 1 ; 
break; 

case 4: 

wO-(l-xl); wl-(xl-x4); w2-(x4-x3); w3-(x3-x2); w4-x2; 

new [ 0 ] - ( wO*cbl [ 15 ] +wl*cbl [11 ] +w2 *cbl [ 8 ] +w3*cbl [2 ] +w4 *cbl [ 0 ] ) /l ; 

new [ 1 ] - ( wO *iribl [ 1 5 1 +wl *mbl ( 11 ] +w2 *mbl [ 8 ] +w3*iribl [ 2 ] +w4 *mbl [ 0 ] ) / 1 ; 

new (2 ] - ( w0*ybl [ 15 ] +wl*ybl [ 11 ] +w2 *ybl [ 8 ] +w3*ybl [2 ] +w4 *ybl [ 0 ] ) 11 ; 

new [33- ( wQ*kbl ( 15 ] +wl*kbl [ 11 ] +w2 *kbl [ 8 ] +w3*kbl [2 ] +w4 *kbl [ 0 1 ) /l ; 

break; 

case 5: 

wO-(l-xl); wl-(xl-x3); w2-(x3-x2); w3-(x2-x4); w4-x4; 
new[0]-(w0*cbl[15]+wl*cbl[ll]+w2*cbl[9]+w3*cbl[4]+w4*cbl[0])/l; 
new [ 1 ] - ( wO *mbl [ 1 5 ] +wl *nibl [ 1 1 ] +w2 *mbl [ 9 ] +w3 *ntol [ 4 ] +w4 *tnbl [ 0 ] ) / 1 ; 
new [ 2 ] - ( wO *ybl [ 1 5 ] +wl *ybl [ 1 1 ] +w2 *ybl [ 9 ] +w3*ybl [ 4 ] +w4 *ybl [ 0 ) ) / 1 ; 
new [3] - (wO*kbl [ 15] +wl*kbl [11] +w2*kbl [9 J +w3*kbl [4 ] +w4*kbl [0] ) /l; 
break; 

case 6: 

wO-(l-xl); wl-(xl-x3); w2-{x3-x4); w3-(x4-x2); w4-x2; 
new[0]-(wO*cbl[15]+wl*cbl[113+w2*cbl[9]+w3*cbl[2]+w4*cbl[03)/l; 
new [ 1 ] - ( wO*mbl (15) +wl*mbl [ 11 3 +w2*mbl [ 9 ] +w3*mbl [2 ] +w4 *mbl [0J)/1; 
new [ 2 3 - ( wO*ybl ( 1 5 ] +wl *ybl [ 11 3 +w2 *ybl [ 9 ] +w3 *ybl [ 2 ] +w4 *ybl [ 0 ] ) /l ; 
new [ 3 ] - ( wO* kbl [ 1 5 ] +wl *kbl [ 1 1 ] +w2 *kbl [ 9 ] +w3 *kbl [ 2 3 +w4 * kbl [ 0 J ) /l ; 
break ; 

case 7: 

w0-(l-x2); wl-(x2-xl); w2-(xl-x3); w3-(x3-x4); w4-x4; 
new[0]-(w0*cblfl5]+wl*cbl[12)+w2*cbl[10)+w3*cbl[4]+w4*cbl[0].)/l; 
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new [ 1 ] - ( wO*mbl 1 1 5 J +wl*mbl (12]+w2 «mbl [ 10 ] +w3 *mbl [ 4 ] +w4 *rabl [0] 
new [2 J - (v0*ybl [ 15 J +wl*ybl [ 12 ] +w2 *ybl [ 10 ] +w3*ybl [ 4 ] +w4 *ybl [ 0 3 ) /l ; 
new[3 J- (w0*kbl [15] +wl*kbl 112] +w2*kbl [10] +w3*kbl [4] +w4*kbl [ 0 J ) 11; 
break; 

case 8: 

w0-(l-x2); wl-(x2-xl); w2-(xl-x4); w3-<x4-x3); w4-x3; 

new ( 0 ] - (wO*cbl [ 15] +wl*cbl [ 12 ] +w2*ebl [ 10 ] +w3*cbl [ 3 ] +w4 *cbl [ 0 ] ) /l ; 

new[l]- (wO*inbl I15]+wl*n»bl [12J +w2*mbl [10] +w3*mbl [3] +w4*inbl [0] ) /I; 

new [2 ] - (wO*ybl [15] +wl*ybl [ 12 ] +w2*ybl [ 10 ] +w3*ybl [ 3 ] +w4 *ybl [ 0 3 ) /l ; 

new [3]- (wO*kbl [15] +wl*kbl [12] +w2*kbl [10] +w3*kbl [3] +w4*kbl [0] ) /l; 

break; 

case 9: 

w0-(l-x2); wl-<x2-x4); w2-(x4-xl); w3-(xl-x3); w4-x3; 
new[0]- (w0*cbl(15]+wl*ebl [12 ] +w2*ebl [ 6) +w3*cbl [3] +w4*cbl [0] ) /It 
new [ 1 ] - (wO*nbl [ 15 J +wl*nbl [ 12 ] +w2*aabl [ 6 3 +w3*rabl [ 3 ] +w4 *nbl [ 0 ] ) /l ; 
new [ 2 ] - (w0*ybl [ 15 ] +wl*ybl [ 12 ] +w2*ybl [ 6 ] +w3*ybl [ 3 ] +w4 *ybl [ 0 ] ) /l; 
new [ 3 ] - (w0*kbl [ 1 5 ] +wl*kbl [ 12 ] +w2*kbl [ 6 ] +w3*kbl [ 3 ] +w4 *kbl [ 0 } ) /l ; 
break; 

case 10: 

w0-(l-x2); wl-(x2-x4); w2-(x4-x3); w3-(x3-xl); w4-xl; 

new [ 0 ] - (wO*ebl [ 15 ] +wl*ebl [ 12 ] +w2 *cbl [ 6 ] +w3*cbl [ 1 ] +w4 *cbl [ 0 ] ) /l ; 

new[l]- (wO*mbl (15]+wl*mbl [12]+w2*robl [6] +w3*robl [l]+w4*nbl [0] ) /l; 

new[2]-{w0*ybl [15]+wl*ybl[l2)+w2*ybl [6]+w3*ybl [l]+w4*ybl [0] ) /It 

new [ 3 ] - (w0*kbl [ 15 ] +wl*kbl [ 12 ] +w2*kbl [ 6 ] +w3*kbl [ 1 ] +w4 * Jebl [ 0] ) /l ; 

break; 

case 11: 

w0-(l-x2); wl-(x2-x3); w2-(x3-xl); w3-(xl-x4); w4-x4; 
new[0]-(w0*cbl[15]+wl*cbl[12]+w2*cbl[7]+w3*cbl[4]+w4*cbl[0])/l; 
new [ 1 ] - ( wO*mbl [15] +wl*mbl [ 12 ) +w2*mbl [ 7 ] +w3*mbl [ 4 ] +w4 *mbl [0 ] ) A; 
new [ 2 ] - (wO*ybl [15] +wl*ybl [12] +w2*ybl [ 7 ] +w3*ybl [ 4 ] +w4*ybl [ 0 ] ) /l; 
new [ 3 ] - { wO*kbl [ 1 5 ] +wl*kbl [ 12 ] +w2 *kbl [ 7 ] +w3 *kbl [ 4 ] +w4 *kbl ( 0 ] ) / 1 ; 
break; 

case 12: 

w0-(l-x2); wl-(x2-x3); w2-(x3-x4); w3-{x4-xl); w4-xl; 

new [ 0 ] - ( w0*cbl [ 1 5 ) +wl*cbl [ 12 ] +w2*cbl [ 7 ] +w3 *cbl [ 1 ] +w4 * cbl [ 0 ] ) / 1 ; 

new [ 1 ] - (w0*mbl [ 1 5 ] +wl*mbl [ 12 ] +w2*rabl [ 7 ] +w3*mbl [ 1 ] +w4 *mbl [ 0 ] ) /l ; 

new [ 21 - (w0*ybl [15] +wl*ybl [ 12 ] +w2*ybl [ 7 J +w3*ybl [ 1 ] +w4 *ybl [ 0 ] ) /l ; 

ne w [ 3 ] - ( w0*kbl [ 1 5 ] +wl *kbl [ 1 2 ] +w2 *kbl [ 7 ] +w3 *kbl [ 1 ] +w4 * kbl [ 0 ] ) / 1 ; 

break; 

case 13: 

wO-(l-x3); wl-(x3-xl); w2-(xl-x2); w3-(x2-x4); w4-x4; 

new [ 0 ] - ( wO*cbl [ 1 5 ] +wl*cbl [ 13 ] +w2*cbl [ 9 ] +w3 *cbl [ 4 ] +w4 *cbl [ 0 ] ) / 1 ; 

new [ 1 ]- (wO*mbl [ 15 ] +wl*n»bl [ 13 ] +w2*mbl [ 9 ] +w3*n»bl [ 4 ] +w4 *nibl [ 0 ] ) /l ; 

new [ 2 ] - (wO*ybl [15] +wl*ybl [ 13 J +w2 *ybl [ 9 ) +w3 *ybl [ 4 ] +w4 *ybl [ 0 ] ) / 1 ; 

new [ 3 ] - (w0*kbl [ 1 5) +wl*kbl [ 13 } +w2 *kbl [ 9 ] +w3*kbl [ 4 ] +w4*kbl [ 0 ] ) /l ; 

break; 

case 14: 

w0-<l-x3); wl-(x3-xl); w2-(xl~x4); w3-(x4-x2); w4-x2; 

new [ 0 ] - (wO*cbl [ 15 ] +wl*cbl [ 13 ] +w2*cbl [ 9 ] +w3*cbl [ 2 ] +w4 *cbl [ 0 ] ) /l ; 

new [ 1 ] - ( wO*inbl [ 1 5 ] +wl*mbl [ 13 ] +w2 *mbl [ 9 ] +w3*mbl [ 2 ] +w4 *mbl [ 0 ] ) / 1 ; 

new [ 2 ] - ( wO*ybl [ 1 5 ] *wl*ybl [ 13 ] +w2 *ybl [ 9 ] +w3*ybl [ 2 ] +w4 *ybl [ 0 ] ) / 1 ; 

new [ 3 ] - ( wO*kbl [ 1 5 ) +wl*kbl [ 13 ] +w2 *kbl [ 9 ] +w3 *kbl [ 2 ] +w4 *kbl [ 0 ] ) / 1 ; 

break; 

case 15: 

w0-(l-x3); wl-(x3-x4); w2-(x4-xl); w3-(xl-x2); w4-x2; 
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new [ 0 ] - ( w0*cbl t 15 ] +wl*cbl ( 13 ] +w2 *cbl t 5 3 +w3*cbl [ 2 ] +w4 *cbl t 0 ] ) / 1 j 
new[lJ-(wO«T«bltl5J+wl*aibl[13]+w2*inblI5]+«3*Bbl [2] +w4*mbl [0] ) /l, 
new [2 ] « ( w0*ybl [15] +vl*ybl [ 13 ] +w2*ybl [ 5 ] +w3 *ybl [ 2 ] +w4 *ybl 1 0 ) ) / 1 
newI3]-(w0*kblfl5J+wl*kblI13]+w2*Jcbl[5J+w3*)cbl [2] +w4*kbl [0] ) /l 
break; 

case 16: 

w(Ml-x3); wl-{x3-x4); w2-(x4-x2); w3-(x2-xl); w4«xl; 

new [ 0] - ( w0*cbl [ 15 ] +wl*cbl [ 13 ] +w2 *cbl [ 5 ] 4«3*cbl [ 1 ] +w4 *cbl [ 0 ] ) /l; 

new [ 1 ] - < w0 *a»bl 1 1 5 ] +wl *ttbl 1 13 ] +w2 *nfcl [ 5 ] +w3 *mbl [ 1 ] +w4 *mbl [ 0 ] ) / 1 ; 

new[2]- <wO*ybl [15]+wl*ybl[13]+w2*ybl [5] +w3*ybl [1] +w4*ybl [0] ) /l; 

new(3]-(wO*kbl[15J+wl*Jcbl[13]+w2*kbl t5]+w3*kbl [l]+w4*Jcbl [0] ) /l; 

break; 

case 17: 

w0-{l-x3); wl-(x3-x2); w2-(x2-xl); w3-(xl-x4); w4-x4; 

new [ 0 ] - ( wO*cbl 1 1 5 J +wl*cbl [ 13] +w2*cbl [ 7 ] 4«3*cbl [ 4 ] +w4 *cbl [ 0 ] ) /l; 

new[l] - (wO*abl [15J+wl*Mbl [13] +w2*nbl [7 J +w3*mbl [4 ] +w4*abl {0] ) /l; 

new [2 ] - (wO*ybl [ 15 ] +wl*ybl [ 13] +w2*ybl [ 7 ) +w3*ybl [4 ] +w4 *ybl [ 0 } ) /l; 

new[3]-(w0*kbiri5]+wl*kbl[13]+w2*kbl[7]+w3*kbl[4]+w4*kbl[0])/l; 

break; 

case 18: 

w0-<l-x3); wl-(x3-x2); w2-(x2-x4); w3-(x4-xl); w4-xl; 
new[0]-(w0*cbl[15]+wl*cbl[13]+w2*ebl[7]+w3*cbl[l]+w4*cbl[0])/l; 
new[l]- (wO*inbl [15]+wl*inbl [13] +w2*iribl [7] +w3*inbl [1] -M^mbl [0] ) /l; 
new [2]- <w0*ybl[15]+wl*ybl [13] +w2*ybl [7]+w3*ybl [l]+w4*ybl [0] ) A; 
new [ 3 ] - ( wO*kbl [15] +wl *kbl [ 13 ] +w2*kbl [ 7 ] +w3*kbl [ 1 ] +w4 *kbl [ 0 ] ) /l ; 
break; 

case 19: 

w0-(l-x4); wl-<x4-xl>; w2-(xl-x2); w3-(x2-x3); w4-x3; 
new[0]-(w0*cbl[15]+wl*cbl[14]+w2*cbl[8]+w3*cbl[3]+w4*cbl[0])/l; 
new [ 1 ] - <wO*mbl [ 15 ] +wl*mbl [ 14 ] +w2*rabl ( 8 ] +w3*mbl [ 3 ] +w4 *rabl [ 0 ] ) /l ; 
new [ 2] - <wO*ybl [ 15 J +wl*ybl [14 ] +w2*ybl [ 8 ] +w3*ybl [3 ] +w4 *ybl [ 0 ] ) /l; 
new[3]-(w0*kbl[15]+wl*kbl[14]+w2*kbl[8]+w3*kbl[3]+w4*kbl[0])/l; 
break; 

case 20: 

wO-(l-x4); wl«(x4-xl); w2-(xl-x3); w3-(x3-x2); w4-x2; 
new [0]- (wO*cbl [ 15] +wl*cbl [14 ] +w2*cbl [8 ] +w3*cbl [2 ] +w4*cbl ( 0 ] ) /l 
ne w [ U- { wO *nibl [ 1 5 ] +wl *ntol [ 1 4 ] +w2 *mbl [ 8 ] +w3 *mbl [ 2 ] +w4 *nibl 10] ) /l 
new [ 2 ] - ( wO *ybl [ 1 5 ] +wl *ybl [ 14 ] +w2 *ybl [ 8 ] +w3 *ybl [ 2 J +w4 *ybl [ 0 ] ) / 1 . 
ne w [ 3 ] - ( wO * kbl [ 1 5 ] + wl *kbl [ 1 4 ] + w2 *kbl [ 8 ] +w3 * kbl [ 2 ] +w4 * kbl [ 0 ] ) / 1 
break ; 

case 21: 

w0-(l-x4); wl-(x4-x3); w2-(x3-xl); w3- (xl-x2) ; w4-x2; 
new[0]-(wO*cblH5]+wl*cbl[14]+w2*cbl[5]+w3*cbl[2]+w4*cbl[0])/l, 
new [ 1 ] - (wO *mbl [ 1 5 ] +wl *mbl [ 1 4 ] +w2 *mbl [ 5 ] +w3 *mbl [ 2 ] +w4 *mbl [ 0 ] ) / 1 , 
ne w [ 2 ] - ( wO *ybl ( 1 5 ] +wl *ybl [ 1 4 ] +w2 *ybl [ 5 ] +w3 *ybl [ 2 ] +w4 *ybl [ 0 ] } / 1 , 
ne w [ 3 ] - ( wO *kbl [ 1 5 3 +wl *kbl [ 1 4 ] +w2 *kbl [ 5 J +w3 * kbl [ 2 ] +w4 *kbl [ 0 3 ) / 1 
break; 

case 22: 

w0-(l-x4); wl-(x4-x3); w2-{x3-x2); w3-(x2-xl); w4-xl; 
new[03-(wO*cbl[153+wl*cbl[14]+w2*cbl[5]+w3*cbl[l]+w4*cbl[OJ)/l; 
new [ 1 ] - (wO*mbl [ 153 +wl*mbl [ 14 ] +w2*mbl [ 5 ] +w3*mbl [ 1 ] +w4 *mbl [ 0 ] ) /l ; 
new [2 3 - (wO*ybl [ 1 5] +wl*ybl [ 14 3 +w2 *ybl [ 5 ) +w3*ybl [ 1 ] +w4 *ybl [ 03 ) /l ; 
new [3]- (w0*kbl[15]+wl*kbl [14]+w2*kbl [5]+w3*kbl [l]+w4*kbl [0] 
break; 

case 23: 
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w0-(l-x4); wl-(x4-x2); w2-<x2-xl); w3-(xl-x3); w4-x3; 

new [ 0 ] - (w0*cbl [ 15] +wl*cbl 1 14 ] +w2*cbl 1 6) +w3*cbl [3 ] +w4 *cbl { 0 ] ) /l; 

new [ 1 ] - (w0*mbl [ 15 ] +wl*mbl [ 14 J +w2*lnbl r 6] +w3*mbl [3 ) +w4 *inbl 1 0 ] ) /l; 

new 1 2 ] - (wO*ybl [ 15 ] +wl*ybl [ 14 ] +w2*ybl [6] +w3*ybl 1 3 J +w4 *ybl [0J)/1; 

new [ 3 J - ( wO*Jebl I15]+wl*kbl [ 14 ] +w2*kbl [6] +w3*fcbl f 3 ] +w4 *kbl [ 0 ] ) /l; 

break; 

case 24: 

w0«(l-x4); wl-(x4-x2); w2-(x2-x3); w3-(x3-xl) ; w4-xl; 

new (01- (w0*cbl 1 15] +wl*cbl [ 14 ] +w2*cbl [ 6 ] +w3*cbl 1 1 ] +w4 *cbl ( 0 ] ) /l; 

new [ 1 ] - (w0*xnbl 1 15] +wl*ntol [ 14 ] +w2*nfcl 1 6] +w3*mbl [1 ] +w4*mbl [ 0 ] ) /l; 

new [ 2 ] - ( wO *ybl [ 1 5 ] +wl *ybl [ 1 4 ] + w2 *ybl [ 6 ] +w3 *ybl [ 1 ] +w4 *ybl [ 0 J ) / 1 ; 

new [ 3 ] - (wO*kbl [ 15 ] +wl*kbl { 14 ] +w2*kbl [ 6 ] +w3*kbl [ 1 ] +w4 *Jcbl [ 0 ] > /l; 

break; 

} 

return; 



1. A method for matching the color display of at least a first display device (12) and at least a second display device 
(14), the first and second display devices having an input, comprising: 

I creating a transform LUT (28) for converting input color values (R,G,B) to output color values (R'.G'.B') by: 

1) producing a first preselected plurality of input color values ((R,G,B) n ), 

2) using said plurality of input values to display a plurality of color patches (32,34) in said first and said 
second display devices (12,14), 

3) obtaining a colorimetric value (Lab) of each of the displayed color patches in each of the displayed 
devices and using said colorimetric values to create a first (40) and a second (42) model for the first and 
the second devices respectively correlating preselected input color values to displayed colorimetric values 
for each of the two devices, 

I I using the transform LUT (28) to transform any plurality of source color image values (R,G, B) to corresponding 
new image values (R',G',B') before inputting said new image values to the input of the second display device 
(14), and 

III displaying said new color image values on said second display device (14), 
characterized in that 

in creating the transformer LUT (28), step I 3) is followed by the following steps: 

4) inputting to said first and second models (40,42) a second preselected plurality of color values (FL^i^G,,^), 
B„ (1) ), 

5) comparing the output (Lab) of the first model to the output (L'a'b 1 ) of the second model to obtain an error 
signal (E) indicative of the difference between the two output signals, 

6) using the error signal (E) to modify the input color values to the second model and again comparing the 
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output of the first model to the output of the second model to obtain a new error signal, 

7) repeating the process of steps 5) and 6) above until the error signal (E) is a minimum, an 

5 8) using the modified color values to create a transform LUT (28) correlating input values (R,G,B) to modified 

values (R'.G'.B'). 

2. The method according to claim 1 , wherein in performing the step II, first testing every one of the plurality of source 
color image values, and 

10 

a) if the source color image values are included in the LUT (28) obtaining from the LUT new image values 
corresponding to the source color image values, and 

b) if the source color image values are not included in the LUT (28) obtaining by interpolation from adjacent 
is source image values in the LUT new image values corresponding to the source image values. 

3. The method according to claim 1 or 2, wherein between steps I 7 and I 8 there is the additional step of optimizing 
the minimum error. 

20 4. The method according to claim 1 or 2, wherein prior to step I 5, the gamut of the output color values of the first 
model (40) is adjusted to be commensurate with the color values gamut of the second model (42). 

5. The method according to claim 2, wherein the source color image values are RGB values and the display devices 
are CRT displays (12,14). 

25 

6. The method according to claim 2, wherein the source color image values are YMCK values and the display devices 
are 4-color printing devices (112,114). 

7. The method according to claim 2, wherein one of the display devices is a CRT display and one a 4-color proofer 
30 and wherein there is included in the color proofer a color converter for converting RGB color values to YMCK color 

values. 

8. The method according to claim 2, wherein one of the display devices is a CRT display (212) and one a 4-color 
printer (21 4) and wherein there is included in the CRT a color converter for converting YMCK color values to RGB 

35 color values. 

9. Apparatus for generating a transform LUT for converting input color values to output color values, comprising: 

1) digital color image values input means (10), 

40 

2) first and second display model means (12,14) for producing each an output colorimetric tristimulus image 
(Lab) values for input color values, each of said model means having an input and an output, the second model 
input connected to the input means (10), 

45 characterized by 

3) an adder device (60) having a first signal input (56) connected to the input means, a second, correction 
error signal input, a control signal input (78), and an output (62) connected to the second model input, for 
outputting modified color values (R'.G'.B 1 ), 

SO 

4) means (66), connected to the output of said first and second model means, for comparing the tristimulus 
values output of the first and second display models and for producing an error signal (E), 

5) means (92) for testing the error signal (E) to determine if said signal is a minimum error signal, and for 
55 outputting a correction eror signal and a control signal to the adder, said testing means (92) connected between 

said means (66) to compare and said adder device (60), and 

6) means (82) also connected to the input means (10), the adder output, and the means (92) for testing, for 
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receiving the adder output and for generating a transform LUT (28) correlating digital input color values to the 
modified color values appearing at the output of the adder (60) when the means (92) for testing determines 
that the error signal (E) is a minimum. 

5 10. The apparatus according to claim 9, further comprising means (82) for storing the generated LUT 

11. The apparatus according to claim 9, further comprising a gamut mapping means (70) for mapping output values 
from the second model (42) into a color gamut commensurate to that of the first model (12), the gamut mapping 
means (70) positioned between the output of the second model and the means to compare (66). 

10 

1 2. The apparatus according to claim 9, further comprising an optimization means (96) in the error testing means (92) 
for determining if the minimum error is an optimum minimum error. 



15 Patentanspruche 

1 . Verfahren zum Abgleichen der Farbanzeige wenigstens einer ersten Anzeigevorrichtung (1 2) und wenigstens einer 
zweiten Anzeigevorrichtung (14), wobei die erste und die zweite Anzeigevorrichtung einen Eingang aufweisen, 
mit folgenden Schritten: 

20 

I Erstellen einer Transformations-LUT (28) zum Transformieren von Eingangsfarbwerten (R, G, B) in Aus- 
gangsfarbwerte (R\ G', B') durch: 

1) Erzeugen einer ersten vorgewahlten Vielzahl von Eingangsfarbwerten ((R,G,B) n ), 

25 

2) Verwenden der mehreren Eingangswerte zum Anzeigen mehrerer Farbflecken (32, 34) auf der ersten 
und der zweiten Anzeigevorrichtung (12, 14), 



3) Erhalten eines kolorimetrischen Wertes (Lab) jedes der angezeigten Farbflecke aufjeder der Anzeige- 
30 vorrichtungen und Verwenden der kolorimetrischen Werte zum Erzeugen eines ersten (40) und eines 

zweiten (42) Modells fur die erste und die zweite Vorrichtung, wobei jeweils vorgewahlte Eingangsfarb- 
werte mit angezeigten kolorimetrischen Werten fur jede der beiden Vorrichtungen korreliert werden, 



II Verwenden der Transformations-LUT (28) zum Umwandeln einer beliebigen Vielzahl von Ursprungsfarb- 
35 bildwerten (R, G, B) in entsprechende neue Bildwerte (R\ G', B 1 ) bevor die neuen Eingangswerte am Eingang 

der zweiten Anzeigevorrichtung (14) eingegeben werden, und 



III Anzeigen der neuen Farbbildwerte auf der zweiten Anzeigevorrichtung (14), 



dadurch gekennzeichnet, daG 

beim Erstellen der Transformations-LUT (28) die folgenden Schritte auf den Schritt I 3) folgen: 



4) Eingeben einer zweiten vorgewahlten Vielzahl von Farbwerten (R n( i), G n(1) , B n(1) ) in das erste und das 
zweite Modell (40, 42), 

5) Vergleichen des Ausgangssignals (Lab) des ersten Modells mit dem Ausgangssignal (L'a'b') des zweiten 
Modells, urn ein Fehlersignal (E) zu erhalten, das die Differenz zwischen den beiden Ausgangssignalen angibt, 



6) Verwenden des Fehlersignals (E) zum Modifizieren der Eingangsfarbwerte fur das zweite Modell und er- 
50 i neutes Vergleichen des Ausgangssignals des ersten Modells mit dem Ausgangssignal des zweiten Modells, 

urn ein neues Fehlersignal zu erhalten, 

7) Wiederholen des Vorgangs nach den genannten Schritten 5) und 6), bis das Fehlersignal (E) einen Minimum 
ist, und 

55 

8) Verwenden der modifizierten Farbwerte zum Erstellen einer Transformations-LUT (28), die Eingangswerte 
(R, G, B) mit modifizierten Werten (R\ G', B 1 ) korreliert. 
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2. Verfahren nach Anspruch 1 , bei dem bei der Durchf uhrung des Schritts 1 1 zunachst jeder der mehreren Ursprungs- 
farbbildwerte getestet wird, und 

a) wenn die Ursprungsfarbbildwerte in der LUT (28) enthalten sind, Erhalten den Ursprungsfarbbildwerten 
s entsprechender neuer Bildwerte aus der LUT, und 

b) wenn die Ursprungsfarbbildwerte nicht in der LUT (28) enthalten sind, Erhalten den Ursprungsfarbbildwerten 
entsprechender neuer Farbwerte durch Interpolieren aus benachbarten Ursprungsfarbbildwerten in der LUT. 

10 3. Verfahren nach Anspruch 1 Oder 2, bei dem zwischen den Schritten 17 und 18 der zusatzliche Schritt des Optimie- 
rens des minimalen Fehlers ausgefuhrt wird. 

4. Verfahren nach Anspruch I Oder 2, bei dem vor dem Schritt 15 die Stufenleiter der Ausgangsfarbwerte des ersten 
Modells (40) derart angepaBt wird, daB sie mit der Stufenleiter des zweiten Modells (42) ubereinstimmt. 

15 

5. Verfahren nach Anspruch 2, bei dem die Ursprungsfarbbildwerte RGB-Werte und die Anzeigevorrichtungen An- 
zeigen (12, 14) mit Kathodenstrahlrohren sind. 

6. Verfahren nach Anspruch 2, bei dem die Ursprungsfarbbildwerte YMCK-Werte und die Anzeigevorrichtungen Vier- 
20 farbendruckvorrichtungen (112, 114) sind. 

7. Verfahren nach Anspruch 2, bei dem eine der Anzeigevorrichtungen eine Anzeigevorrichtung mit Kathodenstrahl- 
rbhre und eine eine Vierfarbenkorrekturdruckvorrichtung ist, und bei dem die Farbenkorrekturdruckvorrichtung 
einen Farbwandler zum Wandeln von RGB-Fa rbwerten in YMCK-Farbwerte aufweist. 

25 

8. Verfahren nach Anspruch 2, bei der eine der Anzeigevorrichtungen eine Anzeige (212) mit Kathodenstrahlrohre 
und eine ein Vierfarbendrucker (214) ist, und bei dem die Anzeige einen Farbwandler zum Wandeln von YMCK- 
Farbwerten in RGB-Farbwerte aufweist. 

30 9. Vorrichtung zum Erstellen einer Transformations-LUT zum Umwandeln von Eingangsfarbwerten in Ausgangsfarb- 
werte, mit: 

1) einer digitalen Farbbildwerteeingabeeinrichtung (10), 

35 2) einer ersten und einer zweiten Anzeigemodelleinrichtung (12, 14), die jeweils kolorimetrische Ausgangs- 

normfarbbildwerte (Lab) fur Eingangsfarbbildwerte erzeugen, wobei jede der Modelleinrichtungen einen Ein- 
gang und einen Ausgang aufweist, wobei der Eingang des zweiten Modells mit der Eingabeeinrichtung (10) 
verbunden ist, 

40 gekennzeichnet durch 

3) eine Addiervorrichtung (60) mit einem mit der Eingabeeinrichtung verbundenen ersten Signaleingang (56), 
einem zweiten Fehlerkorrektursignaleingang, einem Steuersignaleingang (78) und einem mit dem zweiten 
Modelleingang verbundenen Ausgang (62) zum Ausgeben modifizierter Farbwerte (R 1 , G', B'), 

45 

4) eine mit dem Ausgang der ersten und der zweiten Modelleinrichtung verbundenen Einrichtung (66) zum 
Vergleichen des Normfarbwertausgangssignals des ersten und des zweiten Anzeigemodells und zum Erzeu- 
gen eines Fehlersignals (E), 

so 5) eine Einrichtung (92) zum Testen des Fehlersignals (E), urn zu bestimmen, ob das Signal ein minimales 

Fehlersignal ist, und zum Ausgeben eines Fehlerkorrektursignals und eines Steuersignals an den Addierer, 
wobei die Testeinrichtung (92) zwischen die Vergleichseinrichtung (66) und die Addiervorrichtung (60) ge- 
schaltet ist, und 

55 6) eine ebenfalls mit der Eingabeeinrichtung (10), dem Addiererausgang und der Testeinrichtung (92) verbun- 

denen Einrichtung (82) zum Empfangen des Addiererausgangssignals und zum Erzeugen einer Transformat- 
ions-LUT (28), welche die digitalen Eingangsfarbwerte mit den modifizierten Farbwerten korreliert, welche am 
Ausgang des Addierers (60) anliegen, wenn die Testeinrichtung (92) feststellt, daB das Fehlersignal (E) ein 
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Minimum innehat. 

10. Vorrichtung nach Anspruch 9, ferner mit einer Einrichtung (82) zum Speichern der erstellten LUT. 

s 11. Vorrichtung nach Anspruch 9, ferner mit einer Stufenleiteraufzeichnungseinrichtung (70) zum Abbilden von Aus- 
gangswerten des zweiten Modells (42) in eine Farbstufenleiter, die derjenigen des ersten Modells (12) entspricht, 
wobei die Stufenleiterabbildungseinrichtung (70) zwischen dem Ausgang des zweiten Modells und der Vergleichs- 
einrichtung (66) angeordnet ist. 

10 12. Vorrichtung nach Anspruch 9, ferner mit einer Optimierungseinrichtung (96) in der Fehlertesteinrichtung (92) zum 
Bestimmen, ob der minimale Fehler ein optimaler minimaler Fehler ist. 

Revendications 

15 

1. Procede pour mettre en correspondance I'affichage en couleur d'au moins un premier dispositif d'affichage (12) 
et d'au moins un deuxieme dispositif d'affichage (14), les premier et deuxieme dispositifs d'affichage comportant 
une entree, comprenant : 

20 | la creation d'une LUT de transformation (28) pour convertir des valeurs de couleur d'entree (R, G, B) en 

valeurs de couleur de sortie (R', G', B') en : 

1) produisant une premiere plurality preselectionnee de valeurs de couleur d'entree ((R, G, B) n ), 

2) utilisant ladite pluralite de valeurs d'entree pour afficher une pluralite de points de couleur (32, 34) dans 
25 lesdits premier et deuxieme dispositifs d'affichage (12, 14), 

3) obtenant une valeur colorimetrique (Lab) de chacun des points de couleur affiches dans chacun des 
dispositifs d'affichage et en utilisant lesdites valeurs colorimetriques pour creer un premier (40) et un 
deuxieme (42) modeles respectivement pour les premier et deuxieme dispositifs correlant les valeurs de 
couleur d'entree preselectionnees avec les valeurs colorimetriques affichees pour chacun des deux dis- 

30 positifs, 

II I'utilisation de la LUT de transformation (28) pour transformer toute pluralite de valeurs d'image couleur de 
source (R, G, B) en nouvelles valeurs d'image (R\ G', B') correspondantes avant d'entrer lesdites nouvelles 
valeurs d'image dans Pentree du deuxieme dispositif d'affichage (1 4), et 
35 Ml I'affichage desdites nouvelles valeurs d'image couleur sur ledit deuxieme dispositif d'affichage (14), 

caracterise en ce que 

lors de la creation de la LUT de transformation (28), I'etape I 3) est suivie des etapes suivantes : 

40 4) entrerdans lesdits premier et deuxieme modeles (40, 42) une deuxieme pluralite preselectionnee de valeurs 

de couleur (R n( i), G n (i)> B n (i))' 

5) comparer la sortie (Lab) du premier modele avec la sortie (L'a'b') du deuxieme modele pour obtenir un 
signal d'erreur (E) indicateur de la difference entre les deux signaux de sortie, 

6) utiliser le signal d'erreur (E) pour modifier les valeurs de couleur d'entree dans le deuxieme modele, et a 
45 nouveau comparer la sortie du premier modele avec la sortie du deuxieme modele pour obtenir un nouveau 

signal d'erreur, 

(7) repeter le procede des etapes 5) et 6) ci-dessus jusqu'a ce que le signal d'erreur (E) soit un minimum, et 
8) utiliser les valeurs de couleur modifiees pour creer une LUT de transformation (28) correlant les valeurs 
d'entree (R, G, B) avec les valeurs modifiees (R', G', B'). 

50 

2. Procede suivant la revendication 1 , dans lequel, lors de I'execution de I'etape II, on teste en premier lieu chaque 
valeur de la pluralite de valeurs d'image couleur de source, et 

a) si les valeurs d'image couleur de source sont incluses dans la LUT (28), on obtient de la LUT de nouvelles 
55 valeurs d'image correspondant aux valeurs d'image couleur de source, et 

b) si les valeurs d'image couleur de source ne sont pas incluses dans la LUT (28), on obtient par interpolation 
a partir des valeurs d'image de source adjacentes dans la LUT de nouvelles valeurs d'image correspondant 
aux valeurs d'image de source. 
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3. Precede suivant la revendication 1 ou 2, dans lequel, entre les etapes 17 et I 8, il y a une etape supplemental 
d'optimisation de I'erreur minimale. 

4. Procede suivant la revendication 1 ou 2, dans lequel, avant I'etape I 5, la gamme des valeurs de couleur de sortie 
5 du premier modele (40) est ajustee pour etre en accord avec la gamme des valeurs de couleur du deuxieme 

modele (42). 

5. Procede suivant la revendication 2, dans lequel les valeurs d'image couleur de source sont des valeurs RGB et 
les dispositifs d'affichage sont des affichages a CRT (12, 14). 

10 

6. Procede suivant la revendication 2, dans lequel les valeurs d'image couleur de source sont des valeurs YMCK et 
les dispositifs d'affichage sont des dispositifs d'impression en quadrichromie (112, 114). 

7. Procede suivant la revendication 2, dans lequel un des dispositifs d'affichage est un affichage a CRT et un autre 
is un dispositif de tirage d'epreuve en quadrichromie, et dans lequel on a inclus dans le dispositif de tirage d'epreuve 

couleur un convertisseur de couleur pour convertir les valeurs de couleur RGB en valeurs de couleur YMCK. 

8. Procede suivant la revendication 2, dans lequel un des dispositifs d'affichage est un affichage a CRT (212) et un 
autre une imprimante en quadrichromie (214), et dans lequel on a inclus dans le CRT un convertisseur de couleur 

20 pour convertir les valeurs de couleur YMCK en valeurs de couleur RGB. 

9. Appareil pour produire une LUT de transformation pour convertir des valeurs de couleur d'entree en valeurs de 
couleur de sortie, comprenant : 

25 1) un moyen d'entree de valeurs d'image couleur numeriques (10); 

2) des premier et deuxieme moyens de modele d'affichage (12, 14) pour produire chacun des coefficients 
d'image trichromatiques colorimetriques (Lab) pour des valeurs de couleur d'entree, chacun desdits moyens 
de modele comportant une entree et une sortie, la deuxieme entree de modele etant connectee au moyen 
d'entree (10), 

30 

caracterise par 

3) un dispositif additionneur (60) comportant une premiere entree de signal (56) connectee au moyen d'entree, 
une deuxieme entree de signal de correction d'erreur, une entree de signal de commande (78), et une sortie 

35 (62) connectee a la deuxieme entree de modele, pour produire des valeurs de couleur modifiees (R 1 , G', B'); 

4) un moyen (66), connecte a la sortie desdits premier et deuxieme moyens de modele, pour comparer la 
sortie de coefficients trichromatiques des premier et deuxieme modeles d'affichage et pour produire un signal 
d'erreur (E); 

5) un moyen (92) pour tester le signal d'erreur (E) pour determiner si ledit signal est un signal d'erreur de 
40 minimum, et pour produire un signal de correction d'erreur et un signal de commande vers I'additionneur, ledit 

moyen de test (92) etant connecte entre ledit moyen (66) pour comparer et ledit dispositif additionneur (60), et 

6) un moyen (82) egalement connecte au moyen d'entree (10), a la sortie d'additionneur, et au moyen (92) 
pour tester, pour recevoir la sortie d'additionneur et pour produire une LUT de transformation (28) correlant 
des valeurs de couleur d'entree numeriques avec les valeurs de couleur modifiees apparaissant a la sortie 

45 de I'additionneur (60) lorsque le moyen (92) pour tester determine que le signal d'erreur (E) est un minimum. 

10. Appareil suivant la revendication 9, comprenant en outre un moyen (82) pour memoriser la LUT produite. 

11. Appareil suivant la revendication 9, comprenant en outre un moyen de cartographie de gamme (70) pour carto- 
so graphier des valeurs de sortie du deuxieme modele (42) en une gamme de couleurs en accord avec celle du 

premier modele (1 2), le moyen de cartographie de gamme (70) etant positionne entre la sortie du deuxieme modele 
et le moyen pour comparer (66). 

12. Appareil suivant la revendication 9, comprenant en outre un moyen d'optimisation (96) dans le moyen de test 
55 d'erreur (92) pour determiner si I'erreur de minimum est une erreur de minimum optimale. 
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FIG. 7 
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